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CHAPTER I 


GENERAL INTRODUCTION 
1.1 INTRODUCTION 

The evolution and existence of species has been the 
subject of scientific investigation since the days of Darwin. 
Earlier studies were mainly concerned with experimental 
observations and it is only in the beginning of twentieth 
centuary that attempts have been made to predict the evolution 
and existence of species mathematically. The first major attempt 
in this direction is due to Volterra and Lotka which constitute 
the main theme of the deterministic theory of population dynamics 
in theoretical Biology even today. Over the last fifty years, 
many complex models for two or more interacting species have 
been proposed on the basis of Doth a and Yolterra models by 
taking into account the effects of crowding, age structure, time 
delay, functional response, switching etc. (Holling, 1965; 
Rescigno ,1968; 197V; May ,1971; Maynard Smith ,1974; Gomatam, 1974; 
Freedman et al. ,1976; Cushing ,1976; Fred Brauer, 1977; Tansky, 
1978; Juan, Lin et al. 1978 I Goel et al., 1971). 

It may be noted that lotka-Yol terra model focusses on 
population Interactions at a point in space (habitat) ignoring 
movement (migration) which means a perfect mixing of the species 
in a given region., Mathematically, this is equivalent to 
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assuming that the dispersal rates are sufficiently high ahd 
population in the habitat are well mixed. By assuming so, one 
ignores the essential aspects of species response to environmental 
and ecological changes it encounters in the habitat. Thus, 
lotka Volt err a type models di scribe the situations which 
correspond to only laboratory conditions rather than real 
situations arising in natural environment. It may be noted here 
that even in the laboratory spatial variations may be essential 
for the coexistence 0 f the species, (Huffaker, 1958, 1963). 

In the following , there fore ,an account of the literature 
related to migration of the species and its effects on their 
evolution and co-existence is presented. 

1.2 EEFEC1S OE MIGRATION 

In nature, the real habitats could be both homogeneous and 

heterogeneous as regards to their environmental and ecologies^ 
characteristics are concerned. Heterogeneity can arise due to 
topographical and geographical conditions. Seasonal climatic 
changes can also affect the heterogeneity of the habitat. I* 1 
modern times due to rapid industrialization deforestation and 
other man made projects, our environment and ecology have been 
subjected to unrecoverable stresses destroying the natural 
habitats leading to extinction of certain rare species belonging 
to animal and plant kingdom. Due to these and related effects 
the tendency of the species living in the habitat therefore is 
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to migrate to better suited regions of the habitat for their 
survival and existence. In general the movement of the species 
in the habitat may arise due to factors such as crowding, over 
population, predator chasing prey, refuge and fugitive 
strategies. Migration may be caused due to cross wind in case 
of insects, stream of water in case of fishes, climatic 
conditions, flooding, draught, invasions, colonization etc. 

The initial distribution of the population densities or the 
favourable and unfavourable conditions around the habitat may 
also cause the migration of the species even in the homogeneous 
habitats. 


The migration of the species in general can be studied 
by identifying it with convective and dispersive processes in 
both homogeneous and heterogeneous habitats. Since the habitats 
could be both heterogeneous and anistropic it may be remarked 
that the convective and dispersive abilities of the species may 
be different in different directions. Thus, a typical species 
may have convective and dispersive migration in a particular 
direction due to favourable environmental and ecological 
conditions in that direction but it may have no migration in 
other direction due to various hurdles inhibiting their movement. 

The most general model for the evolution of interacting 
species which include interaction, age structure, time delay, 
convective and dispersive migration with space and time dependent 
characteristics could involve equations of various types, for 
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example differential (partial or ordinary), difference-differential, 
integro-diff erential equations, etc. which may include stochastic 
parameters also. These growth equations are already non-linear 
in nature due to interaction of the species and the migration 
terms can make them further nonlinear if density dependent 
dispersal is considered. To make these models compatible, certain 
type of boundary conditions are to be prescribed. These conditions 
depend upon the surroundings of the habitat i.e. whether it is 
surrounded by hostile environment or favourable environment. 

The first successful attempt to study the migration of 
species mathematically is due to Skellam (1951) by identifying 
the migration of species by random dispersal. Since then several 
investigators studied the effects of migration on linear stability 
of interacting species system by considering Vol terra type prey 
predator and competition models, (Kerner, 1959; landahl, 1959; 

Segal and Jackson, 1972; Hadeler ,1974 ; Gopalsamy , 1977; Maynard 
Smith, 1974; Jacob Jorne ,1977-a). 

In particular, Hadeler (1974) has studied the effects 
of dispersion on ' inear stability of interacting species in the 
one dimensional homogeneous finite habitat for prey predator 
model where the dispersion coefficients are equal and constant 
for the two species. It has been shown that effect of such 
dispersion is to stabilize the equilibrium state under non- 
homogeneous boundary conditions giving stable nonuniform spatial 
patterns. Also, the similar analysis has been carried out by 
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G-opalsamy (1977) for competing species in case of homogeneous 
finite and semi-infinite habitats. It is evident from his work 
that dispersive migration has no effect on the otherwise unstable 
equilibrium state in case of semi-infinite habitat while 
migration may stabilize the equilibrium state in the case of 
finite habitat. For unequal but constant dispersion coeffici- 
ents of the two species, the linear stability analysis has also 
been carried out by Mimura and Nishida (1978), Jacob Jorne 
(1977-a) and it has been shown that possibility of dispersive 
instability may arise in such cases. 

As pointed out earlier the habitats are not only 
heterogeneous but anisotropic also, the effects of variable 
dispersion in one and two dimensional habitats would have 
more implications on the stability and co-existence of species. 

In nature there are sufficient evidences of density 
dependent dispersal of species as presented by Wolfenbarger 
et al* (1974), Krebs et al, (1974). Such effect on the 
stability of equilibrium state has been considered by various 
investigators, (Gurney and Nisbet, 1975, 1976; Gurtin and 
MacCamy, 1977) and it has been shown that it further stabilizes 
the equilibrium state. 

The effects of ecological heterogeneity on prey predator 
system has also been investigated by Comins and Blatt (1974) 
by considering the habitat having both favourable and 
unfavourable regions leading to biased dispersion of the 
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species from unfavourable to more favourable regions, giving 
rise a term similar to convective migration in usual reaction 
diffusion equation. They have shown by numerical calculations 
that such biased dispersion of the species has an important 
stabilizing effect. Such biased dispersal can also arise if the 
species are migrating in a particular direction along a gradient 
of some environmental parameter such as wind in case of insects, 
Morris (1963), and water stream in case of fishes. The mathemati- 
cal equation arising out of such biased dispersal has been 
derived by Me Murtrie (1978). 

Recently Me Murtrie (1978) has given a detailed survey of 
literature regarding various aspects of dispersive and 
convective migration for one and two species system. The effects 
of dispersion in a heterogeneous habitat has also been surveyed 
by Levin (1974, 1976). 

At this juncture it may be pointed out that the migration 
of a particular species in a habitat may also be affected by the 
dispersal of other species (cross dispersal), Kerner (1959). 

Rosen (1977-a) has studied such effect on interacting species 
systems with logarithmic dispersion term in the model. 

Keeping this in view, in Chapter II, a generalized model 
has been proposed by considering density dependent self and cross- 
dispersion and its stability has been investigated both linearly 
and non-linearly (See summary, Chapter II, III and VII). 
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1.3 EFFECTS OF PATCHINESS 

As pointed out earlier the real habitats are heterogene o 
and could he inherently patchy due to discontinuous variations in 
topological, ecological and environmental characteristics. Since 
every continuum can be approximated by the finite sum of the 
sequence of constant functions even the habitat with continuously 
varying characteristics can be visualized as a patchy habitat 
with constant but different characteristics in different patches. 
It may, therefore, be very pertinent to ask, can the patchiness 
in- the habitat confer stability to interacting species system 
with convective and dispersive migration terms? 

In general, the behaviour of the interacting species as 
far as their interaction rates and migration rates are concerned 
could be different in different patches. In particular though 
it may happen that in one patch they have dispersive/convective 
ability while in other patch they may not have so due to certain 
unfavourable conditions in that patch. It can also be visualized 
that in the same patch convective migration may exist along one 
direction in the patch but not in other direction due to 
topographical, ecological and other conditions in the case of 
the two dimensional habitats. It can then happen, in some cases, 
as the patches are interconnected, that one patch may become the 
reservoir for the adjoining patch. 
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It can now be remarked, that the equations governing 
the evolution of the species in eacn patch would be of the same 
form as in the case of single habitat but with different 
parameters in different patches. To find their solutions, 
suitable matching conditions have got to be prescribed at the 
various interfaces connecting the patches. These conditions 
normally are the continuity of population densities and fluxes 
at the interfaces. But due to certain topographical and other 
reasons the interfaces or barriers could be such that the 
densities may be continuous but not the fluxes. 

Effects of patchiness in heterogeneous habitat on the 
stability , co-existence and other related aspects of interacting 
species have been studied by several investigators (Cohen, 1970; 
Levins and Culver, 1971; Horn and Mac Arthur, 1972; Vander 
Meer,l973; Levin, 1974; 1976; Chewning, 1975; Maynard Smith, 1974; 
Kawasaki et al., 1979) by considering the movement of species from 
one patch to other as proportional to the difference of densities 
of the respective species in adjoining patches (passive dispers- 
ion). In particular, Levin .(1974) has noted that over all species 
richness is expected to be higher in patchy environment. It 
has been pointed out by Kawasaki et al.,(l979) that heterogeneity 
of environment and nonlinear diffusive motions play an important 
role in stabilizing the populations of competing species. 
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In the above mentioned studies the inherent dispersive 
ability of species in each of the patches which arises due to 
spatial variation causing density gradient within the patch has 
not been taken into account. Iherefore, these studies may be 
useful only in cases where the patches are fairly narrow. 

Keeping this in view, the effects of patchiness on 
linear and nonlinear stability of the equilibrium state of the 
system has been investigated in Chapters IV to VIII by 
considering various combinations of convection and dispersion 
in each of the patches with suitable matching at various interfaces 
[[ See summary • 

1.4 LIAPUNOV STABILITY 

In the previous two sections, we have described the 
literature regarding the effects of dispersal on the stability 
and co-existence of the species in homogeneous and heterogeneous 
(patchy) habitats which arc mainly related to the study of 
linearized version of nonlinear models. I o get the real insight 
of problems, the nonlinear system as a whole must be investigated. 
It is pointed out here that such coupled nonlinear systems cannot 
be solved analytically and one way to study their qualitative 
behaviour is by using Liapunov second method (LaSalle, 1961$ 

Denn, 1975). This method has also been used even to study the 
linear stability of the equilibrium state of interacting systems 
(Gatto et al., 1977). 
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In population dynamics, to study the nonlinear stability 
of the equilibrium state Liapunov second method has been used 
by several investigators (Goh, 1976, 1977, Harrison 1979; 

Jacob J erne et al. *1977, Hastings ,1978, Hsu$ 1978 /Gilpin >1974)* 

In particular, the nonlinear stability of diffusive 
Lotka Volterra system has been carried out by Jacob Jorne et 
al. (1977), Gopalsany et al. (1980) and it has been shown that 
the otherwise stable system remains stable with positive 
dispersion coefficients under zero flux boundary conditions. 

Harrison (1979) has given a Liapunov function which 
generalizes the functions used by Volterra, Goh and Hsu and 
can be used to study the nonlinear stability of variety of 
models even with functional response. Kawasaki (1979) has used 
the Liapunov second method to study the nonlinear stability of 
the prey predator model in the case of patchy habitats with 
passive dispersion. 

Keeping this in view in Chapters VII and VIII the 
linear and nonlinear stability of the interacting system 
in both homogeneous and heterogeneous habitats has been 
discussed by considering dispersive migration of the species 
with density and space dependence. 
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1.5 SUivili'iAIiY 

In the present chapter, the survey of relevent 
literature has been presented so tnat the work done in this 
thesis can be seen in its proper prospective. 

In Chapter II and III, the linear stability of inter- 
acting species with density dependent dispersal has been studied 
while in Chapter VII the corresponding Liapunov stability for 
linear and nonlinear systems has been investigated. The 
Chapters IV - VI are devoted to study the effects of space 
dependent dispersal of species on the linear stability by 
considering patchiness in habitats. Tne corresponding nonline ar 
analysis has been carried out in Chapter VIII. 

In Chapter II, a generalized model for migrating species 
has been derived and it has been snown that in absence of cross- 
dispersion, an equilibrium state which is stable without 
dispersion is always stable with dispersion provided the 
dispersion coefficients of the two species are equal. However, 
when the dispersion coefficients of the two species are 
different, the possibility of self-dispersive instability arises. 
It has been also pointed out that the cross-dispersion of species 
may lead to stability or instability depending upon the nature 
and the magnitude of the cross-dispersive interactions in 
comparison to the self-dispersive interactions. 
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It has been shown that the self convective movement of 
species increases the stability of the equilibrium state and can 
stabilize an otherwise unstable equilibrium state. She effect 
of cross convection (in absence of self dispersion and self 
convection) is to stabilize the equilibrium state in a prey 
predator model with positive cross dispersion coefficients for 
the prey species. 

Finally, it has been noted that if the system is 
stable under homogeneous boundary conditions it remains so under 
nonhomogene ous boundary conditions. 

In Chapter III, the linear stability of two competing 
and migrating species in two dimensional finite and semi-infinite 
habitats is discussed under time dependent reservoir conditions 
and the conditions for stability of an equilibrium state with 
dispersion and convection derived. 

In the case of prey-predator Volterra model an equilibrium 
state is shown to be stable with dispersion and convection in 
all cases. 

In Chapter IV, the effect of space dependent dispersive 
migration on the linear stability of the equilibrium state of 
the system has been discussed by dividing the habitat into 
finitely many patches (subhabitats) where dispersion coefficient 
in each patch is assumed to be constant but different in different 
patches. This is reasonable because variable dispersion 
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coefficient can be approximated by the sum of finitely many unit 
step functions. It has been shown that the dispersion of 
competing species stabilizes an otherv/ise unstable equilibrium 
state in a patchy habitat under reservoir conditions along the 
boundary of the habitat. However, under flux boundary conditions 
the equilibrium state remains unstable. It has also been noted 
that the equilibrium state which is stable in two dimensional 
habitat may not be stable in a linear one dimensional habitat. 

Chapter Y deals with the effects of convection on the 
linear stability of competition and prey predator models in a 
two dimensional patchy habitat with transverse dispersion in 
all the patches. It has been shown for the competition model 
that with zero convective velocity of the species in a particular 
patch, the condition for stability is the same as in the case 
of one dimensional linear habitat. However, with nonzero 
convective velocity in a given patch, the equilibrium state may 
be stable or unstable depending upon the stability or instability 
in the preceeding patches. But if longitudinal convective 
velocity is zero in all the preceeding patches the equilibrium 
state is always stable. 


In case of prey predator model, however, the equilibrium 
state is always stable in all the patches. 


In Chapter VI, the linear stability of migrating species 
in a patchy habitat has been discussed by considering various 
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combinations of convection and dispersion in different patches. 

In the case of one dimensional habitat with dispersion 
only in the last patch, it has been shown for competition model 
that an inherently unstable equilibrium state becomes stable in 
all patches in absence of dispersive migration. However, in the 
case of semi-infinite habitat the migration has no effect on the 
stability of the equilibrium state even in the presence of 
dispersion in the last patch. But when this patch is finite, 
a general condition for stability involving convective velocity 
and dispersion coefficient has been obtained. 

In the case of two dimensional habitat consisting of 
only two patches, the linear stability of interacting system 
has been discussed with dispersion only in the last patch while 
transverse dispersion in both the patches. In the case of 
competition model, it has been shown that the stability of 
equilibrium state in the second patch is dependent on the 
conditions of its stability in the first patch and even if the 
equilibrium state is stable in the first patch, it may not be 
so in the second patch. But if the equilibrium state is 
unstable in the first patch, it would always be unstable in 
the second patch. 

However, in the presence of equal longitudinal convection 
and dispersion of the species in all the patches, it has been 
shown that the effect of longitudinal convection is stabilizing 
in absence of transverse convection. / 
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In Chapters VII, VIII, the linec^ and nonlinear 
stability of equilibr ium state of interacting species system 
in a two dimensional finite habitat has been investigated for 
density and space dependent dispersal respectively by using 
Liapunov second method* It has been shown that conditions 
obtained for the linear stability are both necessary and 
sufficient* further, the condition and region for nonlinear 
stability have also been derived in cases of competition and 
prey predator models. 

finally, it may be remarked here that the method 
developed and used in thesis can provide a new direction for 
studying non-linear partial differential equations with variable 
coefficients. 
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CHAPTER II 


EFFECTS OF CONVECTIVE AND DISPERSIVE INTERACTIONS 
ON THE LINEAR STABILITY OF TWO SPECIES SYSTEM 

2.1 INTRODUCTION 

As pointed out in Chapter I, the classical Vol terra 
model for the evolution of two interacting species ignores 
the effects of migration which may arise due to environmental 
and ecological gradients in the habitat. These effects have 
been studied by taking into account the dispersive and 
convective migration terms in population models and the 
conditions for the linear stability of the equilibrium state 
derived for constant migration characteristics. (Skellam, 1951; 
Landahl, 1959; Segal and Jackson, 1972? Levin, 1974? Hadeler, 
1974; Gopalsamy, 1977; Comins and Blatt, 1974). The 
importance of density dependent dispersal in the case of single 
species model has also been studied (Gurney and Nisbet, 1975) 
which may be selevent in anisotropic habitats. 

It may be pointed out here that the evolution of a 
particular species not only depends upon the self dispersive 
migrations but is also affected by cross (biased or forced) 
dispersion which may arise due to dispersive migrations of 
other species (Kemer, 1959 # Rosen, 1977, a, b) has studied 



stability of the dynamic system by considering cross dispersion 
of the species, using G-omatam (1974) type interactions and 
generalized lick's law involving logarithmic terms. Similar 
ideas have been exploited in chemical kinetics by Jacob Jome 
(1974, 1975). 

Keeping this in view, in this chapter, a generalized 
model for the interaction of two migrating species has been 
derived by considering self as well as cross dispersion effects. 
The linear stability of such a system is investigated by 
using Hurwitz criteria. Effects of convective migration of 
the species on the linear stability of systems have also 
been studied* 

2.2 MATHEMATICAL MODEL 

The equations governing the evolution of two interacting 
and dispersing species system in a heterogeneous habitat has 
been obtained by, Skellam (1951), as 

du. 

= ^i^ U l ,U, 2^ ”* 7 ( J i^ = ,(2*1) 

where u-^Ug denote the densities of the two species and the 
functions f^u^jUg), i = 1,2 discribe the interactions of 
species. The flux vector JL may be given by 
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where 5^ is the constant dispersion coefficient of the i^ 
species. Gurney and Hisbet (1975) have proposed that if the 
transport properties of species are nonlinear, the flux 
may he generalized as 

Ji = -(5 ± + c ijL u.)vu. i * 1,2 (2.3) 

is a constant. The second term may he 
interpreted as the dispersive interaction of the species with 
itself and may be termed as self dispersive interaction. 

The evolution of species may also he affected hy cross 
dispersive interactions (Kerner, 1959). To explain this 
point let us consider the prey predator model, let u^Ug 
denote the population of preys and predators respectively. 

In real situations the tendency of preys would he to keep away 
from the predators and as such the escape velocity of preys 
may he taken as proportional to the dispersive velocity of 
predators (Kerner, 1959). Also the tendency of predators would 

I 

he to get closer to the preys and hence the chase velocity of : 

i 

predators may be assumed to be proportional to the dispersive ' 

velocity of preys. Keeping this in view, the following possibilities 
may arise. ! 

(i) If the preys are escaping away with a velocity which is 

proportional to the dispersive velocity of the predators i.e. J 

vug and are having no encounters with the predators, their 
population is likely to increase.. Therefore a term proportional j 
to vUg, (being negative) should he subtracted from the flux j 


vector Jh 


where c.. 
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Sector J^ of preys [[see equation (2.1) and (2,2)[]. Similarly, 
if the predators are dispersing with the chase velocity which 
is proportional to the dispersive velocity of preys then, 
without encounters with preys, their population may decrease 
as they are wholy dependent on the prey population. In such 
a case, a term proportional to vu^ (being negative) should be 
added in the population flux vector Jg of the predators. 

(ii) While migrating with dispersive velocity proportional 
to v Ug, the predators may have a possibility of encounters 
with preys which is likely to increase the predator population. 
Due to this a term proportional to u^ vug must be subtracted 
from Jg, the flux vector of predators. Since this would 
happen at the cost of prey population, a similar term may be 
added to the prey population flux vector J^. 

(iii) Similarly the preys while migrating with dispersive 
velocity proportional to v u^ may have encounters with predators 
and therefore a term proportional to Ug vu^ should be added 

to j as prey population is likely to decrease. Since this 
kind of interaction may increase the predators population, the 
similar term should be subtracted from the population flux 
vector Jg of predators. 

Taking these points and the interasp'ecies dispersive 
interactions, given by (2.3), into account, the flux vector 
and Jg for preys and predators respectively may be 
written as 



24 


J 1 = -(V 0 llV°12 U 2 ) 'V^lVl 1 VU 2 
S S = - ( V C 22V S 21 U 1 ) , V ( V e l2 U 2 ) 7U 1 


(2.4) 


where 5^, d^, c^» e 12 , e gl , i,j=l ,2 are constant coefficients. 
Similar discussions are also applicable to other types of 
interactions such as competition, 

Thus in general the flux vector of two interacting species 
can he written as follows: 


J 1 = vu 1 - d 1 (u 1 ) vu 2 

Jg = “DgCUfjUg ) VUg - dg(Ug) vu 1 


(2.5) 


It may he noted here that both the flux vectors may depend 
upon the density gradient of both the species and the self 
dispersion coefficients and cross dispersion coefficients 

d^, i=l ,2 may he functions of population densities and can he 
positive and negative depending upon the types of dispersive 
interactions and their magnitudes [see equation (2.4) []• 


Now if the convective migration of both the species are 
also taken into account then it may he noted that due to this 
the fluxes of the two species would increase by 


where U.^ denotes the self convective velocity of ith species 
and U. .i / 3 are convective velocities which might arise due 

1 J 
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to environmental and ecological gradients and biased dispersions, 
(Comins and Blatt, 1974). 

Thus, the flux equation (2.5j may be generalized as follows: 


q = -BiCui.ug) 


TO 2 +B li U l + 72 u 2 

f S = 'VVV 

va 2 -a 2 ( u £ ) 

VU 1 +B 21 U 1 + r, 2Z a Z 


Thus, the equations governing the evolution of two species 
in a two dimensional habitat having environmental and ecological 
gradients can be written from (2.1) and. (2.-6) as 


au, au ± au- au p au p 

~ + U ly ay + V 1 x ~ + V ly ay = i 'l^ u l ,u 2^ + 


W + U - 1 


3x ( D lx ax 


au- ~ au- . 

+ ay ^ D iy ay - ' 1 + ax ^ d ix ax 


au 9 au Q 

<) + <L_ (d £) 

; ay ^ ly ay ; 


au f 


+ u f 


au. 


8u, 


au. 


+ u. 


+ 7, 


3u 


+ V, 


1 


aT" * “2x dx~ * "2y W ‘ v 2x ax” r v 2y ay" ” x 2^ u i*“2' 


= f o(u-,»u-o) + 


a /-n d ~ 2 \ , a__ ^ U 2 \ a_ 

3x ^2x ax ' ay ^ 2y ay ' 6x 


(dr 


au. 


2x 3x 


0 + 


(d 


au-, 


ay v 2 y ay 


(2.7) 

where D. , D. and d. , d. , i=l,2, are the self and cross 
xx xy xx 

dispersive coefficients and U , U. , V. , V. denote the 

xx xy xx xy 

self and cross convective velocities in x and y directions 
respectively . 


The spatially uniform equilibrium state (u^, u g ) of the 
system (2.7) is obtained by 



To study the local stability of the positive equilibrium state 
given by (2.8), the equation (2.7) may be linearized by using 


u. = u. + v. 


i = 1,2, 


(2*9) 


and neglecting the square and higher order terms to get a 


5v n 6v 1 5v 1 dv 9 9v p 2 

W + U lx “ + U ly “ + dir + V ly — = b lk V k 


9 2 v. 


- ^ 

d v-. 


+ D 


lx 3x 2 


+ D 


ly ay 2 


+• d. 


aV 


a 2 

a v r 


131 ax 2 


+ d 


ly ay 2 


( 2 . 10 ) 


ev 9 av 9 av p av 1 av 1 2 

air + U 2x “ + U 2y ” + V 2x “ + V 2y ~ = ^ b 2k v k 


a 2 

d v^ 


.2 2 .2 

a v a v p a v t 

+ f>p„ — 2 ^ + D ?v + d ?x — 2 ~ + d 9 V — 2 

2x dx * 2y dy ~ 2 x Qx c 2y * v 2 


ay 


where 




5f. 

(aTT>- _ . i> k = 

* u ,u 


( 2 . 11 ) 


and v^, i=l,2, are the small perturbations in population of 
the two species. Keeping in view the equation (2.4), it may be 
noted thai the coefficients D.^, ^iy’ d ix* d iy» d=1 >2 are 
functions of equilibrium values u^, u.g, and may be positive or 

5 * 

negative depending upon the type of interactions. 
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The linearized system (2.10) may be associated with the 
following initial conditions 

^(x^y , 0) = T\(x,y) ; i = 1,2 (2.12) 

For the finite habitat of dimensions Lp,Lj x |_0 ,b] , the 
nonhomogene ous reservoir conditions may be prescribed as 

■^(x^t) = ^ il (y) at x = 0 

0 < y < B 

v.j_(x,y,t) = N ig (y) at x = L 

^(xjyjt) = P tl (x) at y = 0 t > 0 (2.13-j 

0 < x < 1 

v jL (x,y,t) = P i2 (x) at y = B 

or the flux conditions (when species are entering the habitat 
from outside) as 



(x,y,t) = M il (y ) 

at 

x = 0 

9x 

(x,y,t) = M ig (y) 

at 

0 < y < B 

x = L 

9v^ 

" dy“ 

(x,y,t) = Q il (y) 

at 

t > 0 (2.13-b) 

y = o 

av ± 

6y 

(x,y,t) = Q ig (y) 

at 

0 < x < 1 

y = b 


In the following section, the local stability of the system 
(2.10) with above initial and homogeneous boundary conditions are 
discussed by ignoring the convective migration of the species. 

The self and cross dispersion coefficients are considered to be 
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constant with respect to the fixed equilibrium state (u^,u g ). 

The coefficients D. , D- are assumed to be positive while 
. »*-*»*. 

d. , d. , i=l,2 can be both positive or negative in all models 

XX J~j 

in contrast to Rosen (197V a,b). 


2.3 EFFECTS OF DISPERSION ON STABILITY UNDER HOMOGENEOUS 
BOUNDARY CONDITIONS 


To study the local stability of the system (2.10) without 
convection, let us write its solution, which satisfies the 
homogeneous reservoir conditions (2.13-a), i.e. N . . = P . . = Oj 

X ^ X 

i,j = 1,2, as follows: 


Vf(x,y ) 


oo OO 

£ £ 
m=l n=l 


xt 


A imn 


sin 


max 

L 


sin 



i=l,2 

(2.14) 


Substituting (2.14) in (2.10) (with convective velocities zero) 
we get the following characteristic polynomial for x 


X - X{t 11 +b g2 -(o 1 +a g )} + b H b 22" b 12 b 2l“ b ll CT 2“ b 22 CT l 


+b 2l P l +b 12 P 2 +Cr l CT 2~ P l P 2 “ 0 


a. = D. 
l ix 


L 

£ 2 
L 

and m,n = 1,2. . . . 


„2 2 2 2 

H- + D iy 


p i 


a , * A S 

a ix “ 2 a iy ^2 


(2.15) 

(2.16-a) 

(2.16-b) 


Similarly, the solution of the system (2,10) (with convective 
velocities zero) under homogeneous flux boundary conditions 



29 


(2.13-b)? i.e. 


M. .= 0 = 

1 D 


Q. 

13 


i ,3 = 1,2? can be written as 


OO 00 

v. = Z Z 
1 m=0 n=0 


At 

e B. 

ion 


cos 


mx 


cos 


nmy 

3 


(2.17) 


which on substitution in equation (2.10) with U. = U. =' 0, 

u v xx 1 y ’ 

V ix = V iy = 0 giYe rise ‘ fco same characteristic polynomial 
as (2.15) with 


to 2 2 J2 

= hx + V ^ 


J 2 2 2 

p = d SL -S— + d ~ 
p i ix l 2 + a iy B 2 


(2.18-a) 

(2.18-b) 


but m,n=0,l,2, . .. 

Using Hurwitz criteria, it may be noted from (2.15) that 
the equilibrium state given by (2.11) is stable, if the following 
conditions are satisfied, 


b ll + b 22 “ (0 1 + °2 } < 0 (2.19-a) 

b ll b 22“ b l2 b 2l~ b ll a 2” b 22 CT l +b 2l P l +b I2 P 2 +<7 i a 2~ P l p 2 > 0 

(2.19-b) 

Bor the case when there is no dispersion the above stability 
criteria simplify to 


b ll +b S 2 < 0 


b U b SS -b 12 b 21 > 0 


( 2 ,20-a) 
( 2 . 20 -.*) 


It is noted from (2,19) and (2.20) that an equilibrium state 
which is stable without dispersion would remain stable provide^ 
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(2*19-b) is satisfied, the possibility of which is real in the 
following cases. 

(i) P 1 P 2 < 0 (2.21) 

(ii) D. x > |a. x i , D iy > |a iy | , i = 1,2 (8.22) 

In a particular case of (ii) i.e. for 


O — CT r=p =s 0 — 0 

1 2 12 


(2.23) 


the otherwise stable state would remain so provided 


b ll +b 22 < b 12 +b 21 


(2.24) 


The following cases are discussed in detail. 
Case 1: Cross dispersion is absent. 


In this case, P i = 0j i = 1,2 and the conditions for 


stability of the equilibrium state simplify to 
b ll +b 22 - (o 1 + Og) < 0 

b ll b 22 _b 12 b 21~ b ll 0 2 _b 22 0 :L +c, l‘'2 > 0 


( 2 . 25-a) 
(2.25-b) 


Considering the case, 


“lx = D 2 x = h > 0 
D ly " B 2y “ D 2 > 0 


( 2 . 26 ) 


i .e, 


CT 1 = 0 2 = a > 0 


it can be noted that if conditions (2.20) are satisfied then 
conditions (2.25) are also valid. Hence, it is remarked that 
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equilibrium state which is stable without dispersion would 
always be stable with self dispersion provided the self 
dispersion coefficients are equal. But when the dispersion 
coefficients are different the condition (2.25-b) may not be 
satisfied even when the conditions (2.20) do indicating the 
possibility of dispersive instability. These results are 
valid both for reservoir and flux boundary conditions. 


Considering the particular case, given by (2.26), under 
reservoir boundary conditions, the criteria for stability (2.25) 
become, 


D>. Dp o 

+ ^ % > b U + b 22 


2(~ + ~) 7t 4 > b? 1 + bf 0 + 2b, Q b 
l 2 b 2 11 


22 


12 21 


(2.27-a) 
(2 .27— b) 


after using the value of o for the least values of m,n as 
unity. The condition (2.27-b) shows that the degree of the 
linear stability of the equilibrium state increases as the 
self dispersion coefficients increase, but it decreases as the 
size of habitat increases. 

It is also noted that the conditions (2.25) can be 
satisfied irrespective of the conditions (2.20) and thus an 
unstable equilibrium state without dispersion may become stable 
with dispersion under reservoir conditions. 

However if an equilibrium state without dispersion is 
unstable it would remain so under flux boundary conditions as 
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the stability conditions in this case are unaffected for at 
least m=0, n=0 £ see equations (2.25) and (2.20) 


In a special case, when b^=0, bg p =0, the corresponding 
characteristic polynomial (2.15) without dispersion (ck=0) 
would have two complex conjugate roots with zero real parts 
for ^12^21 ^ ^ fhe state without dispersion 

is neutrally stable. In such a case, for a^o, i=l,2 
the conditions (2.25) are satisfied and the neutrally stable 
equilibriums state would become stable with dispersion. But if 
bigbgi > 0 the roots of the corresponding characteristic 
equation (2.15) with < I= 0, i=l,2 would have one positive and 
one negative root and the equilibrium state is unstable. In 
such a case the condition (2.25) may be satisfied and an 
unstable state without dispersion would become stable with 
dispersion under reservoir conditions provided 


a il a 21 ~ b 12 b 2l > 0 


where and o Q1 are the values of and <* 0 for m=n=l. 


21 


Gase II : Self dispersion is absent. 

In this case cr i =0, i=l,2 and the conditions for stability can 
be written as 


b ll + b 22 0 (2.28— a) 

b ll b 22“ b 12 b 2l +b 2l P l +b 12 P 2“ P l P 2 > 0 (2.28-b) 

where p. can be both positive or negative. It is noted from 

JL 
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(2.20) and (2.28) that if an equilibrium state is stable without 
dispersion it would remain so with dispersion provided (2.28-b) 
is satisfied for all values of m,n r But it is possible to 
choose m,n such that (2.28-b) is never satisfied for the case 
when and pg are of the same sign and an equilibrium 

state which is stable without dispersion would not be stable with 
cross dispersion. However, if and Pg are of opposite 

signs the condition (2.28-b) may be satisfied and the equilibrium 
state may become stable. 

Further, if the system is unstable without cross dispersion, 
it may become stable under reservoir conditions provided (2.28-b) ■ 
is satisfied and P^, p g are of opposite signs. But with flux 
boundary conditions, an unstable equilibrium would never become 
stable . 


2 .4 PARTICULAR CASES 
2.4.1 PREY- PERLATOR MODEL 


To see the effects of dispersion in the case of prey 
predator model, the following form of functions f^u^jUg), 1=1,2 
are prescribed. 


f l = < a l-Vl - °lV U 1 
f 2 = (" a 2 _b 2 U 2 + °2 U h U 2 


(2.29) 


where a^,b^,c^ are constants and u^,Ug denote the prey and 
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predator densities respectively. She spatially uniform equilibria 
state of the system can be obtained from (2.8) as 


U 1’ U 2 


a 2°l + a l b 2 


a l C 2 


Xi 


: 1 C 2 


b l b 2 * 


C 1 C 2 + b l b 2 


(2.30) 


for a 1 c g - a gb^ > 0. 


She constants b. .j i,j = 1,2 in equation (2.10) can be found 

1 3 

after using (2.11) as 


hi = - b i u i i 


b 12 ° 1 U 1 


b 0 ^ 


C 2 U 2 5 b 22 


* b 2 U 2 


(2.31) 


It may be noted from (2.31), (2.20) and (2.25) that an 
otherwise stable equilibrium state is always stable with self 
dispersion and in absence of cross dispersion under both 
reservoir and flux boundary conditions. 


Since d. , d. , i=l,2 can be both positive and negative 
xx x y 

the following possibilities arise. 


P 1 > 0, Pg < 0 

P 1 < 0 P g > 0 

P 1 > 0 Pg > 0 

P ! < 0 < 0 


(2.32-a) 

(2.32-b) 

(2.32-c) 

(2.32-d) 


Keeping in view (2.32-a) and using (2.31) it is noted that the 
conditions (2.20), (2.19) are always satisfied and the equilibrium 
state without dispersions would always be stable with self as 
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well as cross dispersions. These results are in line with the 
experimental work on prey -predator model by Huffaker (1958), 
Huffaker et al. (1963;. 

It may also be remarked here, by after noting equations 
(2.31), (2.19), (2.20) and (2.32-b,c,d), that the possibility of 
dispersive instability does exist. 


2.4,2 COMPETITION MODEL 


Considering the case of competition model, the interaction 
functions f^(u 1 ,Ug), i=l,2 are given by 

f l = ( a 1 “ b 1 u 1 "' c l a 2^ U 1 

(2.33 ) 

f 2 = ^ a 2 -l3 2 a 2“ C 2 U l^ U 2 


where a. 


.^jb^jC^ are positive interaction coefficients. 


The positive nontrivial equilibrium point (u^,Ug) of 
the system (2.10) in this case is given by 


V U 2 = 


a l b 2 “ °l a 2 
b l b 2 - C l°2 f 


a 2 b l 


°2 a l 


h b 2 - c l°2 


(2.34) 



(2.35-a) 


(2.35-b) 


The constants defined in equation (2.10) for this model after 
using the functions in (2.33) can be obtained from (2,11) as 
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b ll b l U l » b 12 C 1 U 1 
b 21 = “ C 2 U 2 1 b 22 = " b 2 U 2 


(2.36) 


Keeping in view (2.36) find (2.35-a), it is noted from the conditior 
(2.20) that the equilibrium state (2.34) in absence of dispersion 
is stable. She stability of this system is further enhanced in 
presence of self dispersion and without cross dispersion £ see 
equation (2.25) 3 • However under the condition (2.35-b) the 
equilibrium state (2.34 ) is not stable as the condition (2.20-b) 
is violated. In such a case the equilibrium state may become 
stable with self dispersion provided the condition (2.25-b) i-s 
satisfied • 


Considering also the effects of cross dispersion and 
noting the signs of d lx> & gx , d ly , d gy for this model, it is 
seen that either and P g are of opposite signs (one 

competitor is superior than the other and having gains due to 
cross dispersive mobility) or both are negative. Then, in the 
former case the stable equilibrium state in the case of nonzero 


intraspecies interactions and without dispersion may remain 
stable provided (2.19-b) is satisfied. But in the later case 
it would become unstable as (2.19-b) is never satisfied for all 
values of m and n. However, if the equilibrium state is unstable 
without dispersion it may become stable with dispersion provided 
(2.19-b) satisfied in the case of reservoir conditions. But the 
unstable equilibrium state remains unstable for p^, p g with same 
signs under, flux boundary conditions. 
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2.5 EFFECTS OF CONVECTION ON STABILITY UNDER HOMOGENEOUS 

BOUNDARY CONDITIONS 

In real situations the two interacting species may have 
different convective and dispersive utilities and it may be 
reasonable to postulate that the species having more convective 
velocity can disperse faster. It can then happen that the 
ratios of convective and dispersive velocities of the two species 
may be equal. 

To visualize the difference between the effects of this 
relevent case and the case of equal migration characteristics, 
on the linear ' stability of the equilibrium state, the following 
cases are studied. 

i. Equal convective velocities and dispersion 
coefficients of the two species. 

ii. Unequal convective velocities and dispersion 
coefficients with their ratios being equal. 

To avoid complexities the cases of self and cross migration 
are dealt with separately. 

2.5.1 EFFECTS OF SELF CONVECTION AND DISPERSION 

To see the effects of convection on the stability of the 
equilibrium state, system (2.10) in absence of cross dispersion 
and cross convection can be written as follows 
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9v 1 9v 1 5 t 1 

ar + u ix air + u iy ay 1 


Sv p 0v 9 9v 9 

ar^air* U 2 y w- 


? 2 

2 3 9 v* 

2 b lV v. + D- — jr 1 + D — s- 

3e=i llc k lx ax 2 ly ay^ 


2 2 
a v 9 a Vg 

+ D sx 7^ + D sy “2" 


(2,37) 


2 

2 

k=l 


ay 


where h k=l , 2 is defined as in (2,11). 

Considering the case, 


Ik = Ik = Ik , Ik = Ik = IT y 
lx 2x 1 ly 2y 2 


D lx D 2 x D 1 * D ly 1) 2y ' D 2 ’ 


(2.38) 


the solution of the system (2.37) with (2.38) satisfying the 
homogeneous reservoir conditions (2.13-a) (i.e. N . . = P. - = 0; 

J- J J- J 

i,j - 1»2) can he written as 


QO OO 


v i “ V ex P C 2B^ x + S~ y - xt 3 sin(2|Z.) 

i = 1,2 (2.39) 

Substitution of (2.39) in the system (2.37) give rise to the 
following characteristic polynomial in X 


X 2 - X{b 11 +bgg-2a>+ b-, 1 b o0 -b 1 0 b 01 -(b 1 1 +b 99 ) a + a = 0 (2.40) 


11 22 12 21 v 11 22 


where 


b. - = b 


U? U 2 

*"( AT? ~ + 317“ ) 


11 “ 11 V, 4D- 4D 


b 22 “ b 


U 1 

22 “ ( 4D7 + 4D ^ 


(2.41) 





which gives the conditions for stability of the equilibrium 
as follows: 


b ll +b 22" 20 < 0 


b U b 22~ b 12 b 21~( b ll +b 22^ a + a > 0 


(2.42) 


It can be seen that for TJ^Ug^Q condition is 

equivalent to the conditions (2.25) with a ±- c! 2 =a * can 
therefore be noted from (2.42) that if the system is stable 
without convection it will remain stable with convection also 
and the degree of stability increases as the magnitude of 
convective velocity increases. It is further noted from (2.42) 
that if the system is unstable without convection it may become 
stable with convection provided (2.42) is satisfied. 

Considering the second case, when the convection and 
dispersion of the two species are not equal but the ratios of 
dispersion coefficient to the convective velocity for the two 
species are equal i.e. 


Zlx _ ^2x _ ^2£ 

D lx D 2x ’ D ly D 2y 


(2.43) 


the solution of the system (2.37) satisfying homogeneous reservoir 
ponditions (2.13-a)can be written as 


oo oo 


U. IT 

v ± (x,y,t) = 2 £ Af exp x + 

1 m=l n=l 111111 ^ix iy 


y +Xt] 


sin(S^S) s ixi(~^2L) 
i = 1,2. 


(2.44) 
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which on substitutiong in the system (2.3?) give rise to the 
following characteristic polynomial in 


X " X{b ll +]:, 22 ^ 0 1 +CT 2^ + b ll 1:3 22“ b 12 l3 2l“ < ' ]:) ll Cr 2 +b 22 CT l^ + ^ 0 


12 


where 


b. . 
H 


ii 


IT 2 

IX 

4D . 
ix 


u 2 

-i2L 

4D- 

iy 


(2.45) 


(2.46) 


The conditions for stability in this case can be derived 
as follows: 


b ll + b 22 ~(°1 + °2 ^ < 0 
b ll b 22"" b 12 b 21"*^ b ll a 2 +b 22 CT l^ + a l 0 2 > 0 


(2.4?) 


Since the condition (2.4?) is of similar form as that 
of (2.42), it can be remarked here that unequal convection 
also stabilizes the equilibrium state. Thus, it may be 
remarked that the difference between tho two above mentioned 
cases are only quantitative in nature and qualitatively convection 
(in both the cases of equal and unequal convection) has stabilizing 
effect. 


2.5.2 EPPECTS OP GROSS CONVECTION AND DISPERSION 

Let us write the system (2.10) in absence of self 
convection and self dispersion as 
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dv. dv p 3v 2 a 2 v p 8 2 v 

a T" + v lx ~ + v ly ay“ = + d lx “^2“ + d ly 


(2.48) 


av„ 


oV 


av 0 av- 

a”t*" + V 2x dx" + V 2y ay” “ ^^ b 2k v ]£ + d 2x ^2 

Considering a particular case of (45), 


a v- 


+ d 


2y dy 2 


1! 

$ 

rtf 

li 

w 

C\2 

'■d 

d l 

d iy 

= + d 2y _ 

d 2 

V- = 
lx 

V 2x = 

V 1 

Y iy 

- + v sy = 

Y 2 


(2.49) 


the solution of the system (2.48) with (2.49) under reservoir 
conditions can be written as 


v. 

l 


2 

m=l 


2 

n=l 


\mn e ^ 


( !l 

'2d, 


x + 


y + At) sin(2j££) sin 


(Sg*) 


(2.50) 

Substitution of (2.50) in (2.48) with (2,49) give rise to the 
following characteristic polynomial in X 


X 2 -x(b 11 +b g2 )+b 11 b 22 -b l2 b 21 +p( b 12 +k 21 ) + p 2 - 0 

where p p 

Vj 1 % 

b 12 = b 12 “ 4d” “ 4d^' 

V 2 V 2 

b 2l = b 2l ~ 4d^ “ 41^ 


(2.51) 


(2.52) 


Thus the criteria for stability of the equilibrium state in this 
case can be written as 
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b ll +b 22 < 0 

b ll b 22~ b 12 b 21 + p ^ b 12 +b 2l' > ” p > 0 


(2.53) 


It can be noted from (2.53) that the equilibrium state 
which is unstable without cross convection remains so with 
cross convection also. 

Considering the next interesting case when the convection 
and dispersion of the two species are not equal but the following 
relations between cross dispersion and convection hold 


v Y 

lx = 2x 

di d„ 

lx 2x 


V V 
1.Y _ _2Z 

d ly d 2y 


(2.54) 


Then,' the solution of the system (2.48) under resevoir 
conditions can be assumed as 


v i (x,y,t) 


m= 


OO OO V. 

S 2 B imn exp x + W" y + 

i=l n=l imn ^ a ix ^ a iy 


sin(-~y sin(2gZ.) 


(2.56) 


which on substitution in (2.45) gives the following characteristic 
equation for X 


X 2 -X(b 11 +b 22 )+b 11 b 22 -b l2 b 2l +(P 2 b* 2 +P 1 b| 1 )-P :L P 2 - 0 


(2.57) 


where 
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b* _ b 


. liL. 

12 “ d 12 

~ 4d lx 

4 S 

b* — b^,-i 

Vo 

v 2x 

Vo 

2y 

2l 21 

4d 2x 

4d. 

a 2y 


(2.58) 


and the conditions for stability in this case are obtained as 


b ll +b 22 < 0 

b + ( P £ b l2 ^ *“ P -i P 9 ^ ® 


(2.59) 


1 21 ' 


12 


Keeping in view of (2.53), it can be noted from (2.59) 
that when the cross dispersion coefficients are of opposite 
signs as in the case of preypredator model (Kerner, 1959), the 
equilibrium state which is unstable with equal cross convection 
and dispersion may become stable with unequal cross convection 
and dispersion. 


2.6 EFFECTS OF NONHOMOGENEOUS BOUNDARY CONDITIONS 

To see the effects of nonhomogeneous boundary conditions 
(2.13) on the local stability of the system (2.10), we prove the 
following theorem. 

Theorem s If the linearized system (2.10) with homogeneous 
boundary conditions is stable it would remain so even under 
time independent nonhomogeneous boundary conditions (2.13), 

Proof : Let v^(x,y,t), i = 1,2 be the solution of the system 
(2.10) with the initial condition (2,12) and the boundary 
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conditions (2.13-a) or (2.13-b). let v^(x,y) be the steady 
state solution of the corresponding system (2.10) satisfying 
nonhomogeneous boundary conditions (2 a 13 a) or (2.13-b). Then, 
it can be noted from (2.10) and its corresponding steady 
state that 

z ± (x,y,t) = v^(x,y ,t ) - ^(xjy), i = 1,2 (2.60) 

would also satisfy the system (2,10) with the following initial 
conditions, 

z^x^O) = v i (x,y,0) - ^(xty), i = 1,2. (2.61) 

and the homogeneous boundary conditions corresponding to (2.13). 

Since the forms of z^(x,y,t), i = 1,2 would be the 
same as that of v^(x,y,t), it is noted from equation (2.60) 
that the conditions for stability under both homogeneous and 
time independent non-homogeneous boundary conditions would 
remain the same. 

2.7 CONCLUSIONS 

The evolution and local stability of two interacting 
species in a heterogeneous environment have been investigated 
by taking into account the effects of self as well as. cross 
dispersions and convections under reservoir and flux conditions 
prescribed at the boundary of the finite two dimensional 
habitat . 
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Effects of dispersion (in absence of convection) on the 
stability of the equilibrium state have been summarised in 
fable 1* 

Effects of self convection (in absence of cross dispersion 
and convection) have also been investigated and it has been shown 
that it generally stabilizes the system under both type of 
boundary conditions. In particular, it has been noted that a 
stable equilibrium state without self convection becomes more 
stable while an unstable equilibrium state may become stable 
with convection. 

Effect of cross convection (in the absence of self 
dispersion and self convection) on the prey -predator model is 
to stabilize the equilibrium state. 

Further, it is noted that for equal dispersion coeffi- 
cients and convective velocities the equilibrium state remains 
unstable jjse° equation (2.53)3] • 

A particular case of unequal convection and dispersion, 
when the ratios of dispersion coefficients to the convective 
velocities for the two species are equal, has been discussed and 
it has been shown that the equilibrium state which is unstable 
with equal self and cross dispersion may become stable in this case. 

It has also been shown that if the system is stable with 
homogeneous boundary conditions it remains so under time 
independent non- homogeneous boundary conditions. 
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CHAPTER III 


STABILITY OB T\*0 INTERACTING AmD MIGRATING 
SPECIES SYSTEM IN TWO DIMENSIONAL HABITATS 

3.1 INTRODUCTION 

In the previous chapter, the effects of self and 
cross dispersion and convection on the linear stability of 
two interacting species system nave been discussed in a two 
dimensional finite habitat and the criteria for stability 
of the equilibrium state derived under homogeneous and non- 
homogeneous constant boundary conditions. It can be noted, 
however, that neither the density distributions of the species 
in the habitat have been found nor the effects of time dependent 
boundary conditions have been investigated. 

In this chapter, therefore, tne linear stability of 
two interacting and migrating species system without cross 
convection and dispersion is carried out in finite and semi- 
infinite habitats under time dependent nonhomogene ous boundary 
conditions. The steady and unsteady state distributions of 
the species are obtained. A special case when the dispersion 
coefficient in a particular direction is zero has also been 
investigated. 
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3.2 BASIC EQUATIONS 


Tne system governing the evolution cf two interacting 
and migrating species in a two dimensional anisotropic habitat 
(0 < x < L, y>0) can be written as 


du. 


8u. 


rr + D i st 


au. ~ 

+ Y i aV” f i< u l >U P + I 


0X 

a_ 

ay 


OU. -J 

■ D ix^ u l ,u 2^ 0x~ J 
au. -j 

_ D iy^ u l ,u 2' ) ay” J ^ 3,1 - ) 


where u^, i = 1,2. denotes the density of ith species at any 
point (x,y ) in the habitat at an instant tj , V. , the 
convective velocity components and D^ x , , the dispersal 

coefficients in x (longitudinal) and y (transverse) directions 
respectively which may be, in general, functions of u^,Ug. The 
function ^j_( u i » u 2 1 = 1»2 denotes the interaction term for 

ith species. 


The nontrivial positive equilibrium point of the system 
(3.1) is given by 


f i (u 1 ,Ug) = 0; i = 1,2 (3.2) 

Keeping in view of (3.2) and writing, 

u i (x,y,t) = u i + v j _ (x,y,t) (3.3) 

the system (3.1) can be linearized about the equilibrium state 


as follows : 
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8v. 


dv. 


av i 


2 

8 v. 


FT + D i 53T + v i aT * ^ b ik T k + D il u i’V 

a 1 . 


where 


b - !£i 
ij " au. 


+ J,yl,^ 1 »S 2 ) —g 


i, 8 = 1»2 ♦ 


(3.4) 


U 1’ U 2 


It is noted here that the dispersion coefficients in (3.4) are 
functions of (u^,Ug) and are constant for fixed equilibrium 
state . 


Ihe following initial and boundary conditions for the 
system (3.4) can be prescribed, 


v^(x,y,0) = 0; i = 1,2 ; 0 < x < 1, 0<y<B (3.5) 

^(O.yjt) = 0; i = 1,2 ; 0 < y < B, t > 0 (3.6) 

■v^tljyjt) =' 0; i = 1,2 ; 0<y<B, t>0 (3.7) 

v 1 (x,0,t) = N(t) 

0 < x < 1; t > 0 (3.8) 

v 2 (x,0,t) = P(t) 
such that 

Lt N(t) = K rt , Lt P(t) = P n (3.9) 

t - OO 0 t - CO 0 

v^XjBjt) = 0; i = 1,2; 0<x<L;t>0 (3.10-a) 


for finite habitat 


l i J . f\ 

* •^ r % i ha, l 


m. 
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or 


Lt v-(x,y ,t) = 0; i = 1,2 ; 0 < x < Lj t > 0 (3.10-b) 

y -* oo x 


for the semi infinite habitat in y direction. 


As in the previous chapter, the linear stability of 
the system (3.4) is investigated for the competition and prey 
predator models in the following two cases : 

1. When the dispersion coefficients and convective velocities 
of the two species are equal. 

2. When the ratios of dispersion coefficients to the convective 
velocities for the two species are equal. 

3.3 COMPETI HON MODEL 

In the case of competing species the interaction 
functions f^, fg may be prescribed as 


f — u ( a^ 


C ]_ u 2 ^ 


(3.11) 


fg ™“ Ug (a.g *“ c g i u ^ ) 

and the linearized system for the evolution of two species 
would be given by (3,4) with 


b ii - 0 » i - 

b 12 = = "■ c i a 2/ c 2 0 ( 3 *12 ) 


i 


^21 *""* "" c 2 a i/ ^ ® 
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3,3.1 EQUAL COMVE'GTIOij AND DISPERSION 


Considering rue case 


11 ix 


Jj . 

xy 


rJ i 


v i 


Di 

-L. 

D 2 

U 

V 


1,2 


(3,13) 


and using the nonsingular matrix transformation ^ 0, Dg ^ 0 
for u i f °> u g f 0) 


X = 


I u V ,U 2 V 2 \+"l 

eXp | 257 X + 2D 0 y “ ^4D“ + 4Dg 


[ 


AX. 


(3.14) 


X 1 = exp 


2 

i- "^1 


c£: + V *- <gBT* + A_1 


where , 



1 v 1 (x,y,t)“ 


w 1 (x,y,t)“ 

X = 


o y = ' 

9 ^2. ! 



L v 2 ( x,y,t ^- 


L w 2 ^ x,y,t ^- 


A = 


r cosh (st) sinh (st) 

-sinh (st) -1/r cosh (st) 


r = fI7g ; s = fig 

the system (3.4) with (3.12) reduces to the following form, 
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aw i 

at 




+ d. 



i = 1,2 


(3.15) 


The initial and boundary conditions (5.5) - (3.10) are 
transformed as follows : 


w i (x,y,0) = 0 
w i (0,y,t) = 0 
w i (L,y,t) = 0 


(3.16) 

(3.17) 

(3.18) 


W 


iCxjOjt) = exp £- 


Ux 

W. 


-+ 


(U 2 


)tj S(t)/r cosh (st) 


Wg(x,0,t ) = exp 



+ P(t) sinh 



tt2 

( 4D^ + 


-ff(t) sinh (st) 


w^(x,B,t) = 0 


- P(t ) r cosh 



(3.19) 


or (3.20) 

It w i (x,y,t) = 0 ; i = 1,2 


(i) Finite habitat : 

Solving (3.15) with conditions (3.16) - (3.20) in the 
case of a finite habitat (0 < x < 1, 0 < y < B) by using double 
Fourier finite sine transform with respect to x and y and then 

inverting the resulting solution with the help of (3.14), we get 
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4Dg7l 

v,(x,y,t) = ■ exp 

1 LB 2 


U V 

l 2D 1 x + 2B f y 


(3.21) 


oo oo 


. £ 


m=l n=l 


n G- - (t ) e(m) sin (Sp) sin (^p)s i=l,2 


"E 


where , 


G-i(t) = / V il (I) exp ^(n^njaQdO! 

2 2 D.m 2 ^ 2 D p n 2 7t 2 

I , 1 (m,n) = 2jp + ij- + - — r — + — g — 
^1 ^2 L 2 B 2 


(3.22) 


(3.23) 


0(m) = / exp ( T ° 


g-Q-^ x) sin (SH)dx 


o 


V 11 (T) = N( t) cosh (si) - P(t)r sinh (si) 
Vgl(T) = P(t) cosh (si) - N(t)/r sinh (si) 


(3.24) 


(3.25) 


Keeping in view of (3.9) end (3.25) it can he noted 
that the finite limit of v^(x,y,t), i = 1,2 i.e. v?(x,y), exists 
provided 


B^(m,n) - s > 0 


(3.26) 


and is given hy 


v l (x>y} = £^2 ^ [5§£ x + 25£ y ] 

OO OO kM *— | 

Z £ P 1X [2' 1 (m,n)N 0 -E 0 f J sin(^)sin(£§£) (3.27-a) 
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v«(x,y) = exp 
^ LB^ 


U 

l.2D 1 


x + 


V 

2s; 


y 


] 


CO QO 

2 2 B 

m=l n=l 


U jp 1 (m,n)P 0 -S 0 g 


sin (^) sin (£f4 

(3.27-b) 


where 


B 


11 


B?(m,n) - s^ 


] 


-1 


It can he verified that (3.27) is also the steady state solution 
of the corresponding system leading to the asymptotic stability 
of the equilibrium state under the condition 


o g 

-y-2 y2 D^TC 

+ + — + ir 


s > 0 


(3.28) 


obtained from (3.26) for m = n = 1. 


It is pointed out here that the stability condition 
(3.28) is the same for both the constant and time dependent 
boundary conditions, provided (3.9) is satisfied. 

Burther, even if the condition (3.28) is not satisfied 
an equilibrium state would be asymptotically stable under the 
condition 

N » Pr (3.29) 


to the corresponding steady state distribution given by 


T i (x ’ y) = ~2 exp [Sq + 

l 2 N 0 B 1;L (m,n)sin(^) sin (^) (3.30~a) 
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» 





] 


CO CO 

2 2 
m=l n=l 


F u (m,n) sin (2p) sin (^-) 


•nity 


(3.30-b ) 


where 



P 1 (m,n) 


+ 



( ii ) Semi-. infinite habitat ; 

In the case of a semi- infinite habitat (0 < x < 1, 

0 < y < °°), the solution of the system (3.4) with initial and 
boundary conditions (3.5) - (3.10) can be obtained by taking 
first Fourier finite sine transform with respect to x and later 
Fourier sine transform with respect to y of the system (3.15) and 
then inverting with the help of (3.14), we get 

v j_(x,y s t) = ~ exp . 


where 


2 0 (m) f 


t -3/2 
T 


o 


'7 il (l) exp 



df 


(3.31) 


F 2 (m) 




2 P 

D 1 m k 
l 2 


(3.32) 


As t -* ®, it can be noted from (3.31) that the limit of 
v i(x,y,t) would exist provided 


- s > 0 


(3.33) 
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and is given by 


v®(x,y) = ^ exp 


H 

U 3D 1 


x + SET y 


CO 

y, 0(»>)| [(w> 

ilL=X u ~ 


exp (- ~~) + (H 0 +P 0 r)exp (- sin (222) 

YP 2 0 0 yd 2 -J L 


(3.34) 


v|(x,y) - | exp x + |^yj e(m) | £(P 0 -V r) 


P 1? y P 9P y -j 

exp (- - = — ) + (P 0 +N 0 /r) exp (- Jsin ( 


fV, 


fv< 


■mvL } 


where 


F 12 “ 


-p> _2 2 

I^m ti 


U 2 V 2 

®T ^ l 2 


- s 


2 2 


- 2 U 2 V s D l m It 

J 22 = £57 + fE7 + TV 


+ s 


(3.35) 


(3.36) 


which are the same as the steady state solution showing the 
asymptotic stability of the equilibrium state under the condition 

T j 2 v 2 D ± % 2 

Wl + 3I£ + “ s > 0 (3 * 37) 


Again, in this case also the equilibrium state may 
become asymptotically stable under the condition (3.29) 
irrespective of the condition (3.37), to the steady state 
solution given by 
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= !■ exp 



00 r f 9py-i 

2 e(m) H exp - —z— 
m = 1 Y" Dg ™ 


- n r mnx ) 

iD >i-Xi ) 



which is the same as (3.34) with (3.29) assuring the 
of the steady-state, 

3.3.2 UNEQUAL COKVECilOh AND DISPERSION 
Now consider the case 




V 1 _ V 2 
D ly D 2 y 


which seems to be more realistic 


Using the following transf ormation 


v i = w i 


U. V. 

exp ( 5 * 5 ^- x + y); i = 1,2 

ix iy 


the system (3.4) reduces to 
9w- 


a 2 w. 


at 


ew i - fw 2 + D i x 


+ D 


9 2 w^ 

ly 2 d 


2 2 

3w 2 , 9 Wp 9~Wg 

rr = _sw i - hw 2 + d 2x + D 2y 


(3.38) 


uniqueness 


(3.39 ) 


(3.40) 


(3.41) 
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where 




h = 




(i) Finite habitat 


(3.42) 


ihe solution of the system (3.41; with initial and 
boundary conditions (3.5) - (3.10) under the transformation 
(3.40) can be obtained by using Laplace and double Fourier sine 
transforms with respect to t and x,y respectively from which the 
solution of the original system (3.4) can be written as follows : 




U 1 V 1 ua & |— 

v 1 (x,y,t) = exp (gjji- r + y) £ £ / A ex p(- Pl Ti 

LB^ ^lx ^ly m=l n=l o L lnm lDm 


CO CO “t I 


+ B lmn exp (-P 2mn 2) df sin (2JE) sin (2g£) 


(3.43) 


v 2 (x,y ’ t) = Ttf. exp ( 5UV x + SET" y) E E -T l A Smn exp(-P lmjl I) 

Jj±> cX &V ui~ jl 0 u 


U 2 V 2 


CO C50 t 


r. 


2y iu— X Cl— 1 0 


where 

A 


Iron 


B 


lmn 


+ Bgmn ^(-Pgum 1 '] 41 sln s “ 


B ly ( Plnm ~ H mn> + f B gy f«(m,t-T) 

^lmn ” ^2inn 

- W - f P 2y 

p lmn “ ^2mn 


(3.44-a) 
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B gy (~ p Iran ~ \n,n^ + g ^pj N*(mit-I) 

^Imn “ ^2mn 


^2y^ B mn " Pgan^ Z g D 1.V 

^imn “ ^2mn 


(3 .44-b ) 


$ Imn ’ ^2mn 


-*[> 


mn + V) ± V (V- ii I nn> 2+4 f « 


(3 *45) 


■p p i T) m TU « t\ XI % 

E mn D lx l 2~ D ly ^2 


m2 2 „2 2 

H = h + V„ 2-2- + d EJL. 
mn 2x l 2 2y b 2 


(3.46) 


W*(m,t) = N(t) / exp (- x) sin (2M) dx 

o ^lx h 


P*(m,t) = P(t) / exp (- x) sin (22E) dx 


(3.47) 


It can be noted from (3.43) that as t — the limit of 
v^(x,y,t)j i = 1,2 exist, provided 


tt 2 „2 2 ^2 tt 2 ,,2 2 -p. 2 

f u i , V 1 , ^ly 71 W U 2 V 2 1> 2x 1t ^y 71 ) _2 . 0 

(rrs — + IS — + — + — 9 — ^4D + 4D + 9 — + 9 '~ s > u 

4JJ lx 4 ly L 2 B 2 ^2x ^2y L 2 B 2 


(3.48). 


and is given by 


.8/ „ \ 471 f 1 


v,(x,y) 


exp (■ 


V. go o° f-»A^ 

, 1 \ m V „ r lmn . 2mn 

X + or) y / 2 ^ ^ n n' 

^ly m=l n=l L-Pimn P2mn- 


sin (2|EE.) sin (^^pO (3.49-a) 
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4y__ 

LB^ 


exp ( 


U, 

2D 


2x 


x + 



oo oo 

S I 
m=l n=l 


A 2mn S 2m n 
^ Imn -P 2nm 


sin (H^L) sin (£|Z) (3.49-b) 

where I imn » B inm are the same as A ijn , B inm respectively with 

N(t), P(t) replaced by H Q , P , It is pointed out here that 
(3.49) is the same as the steady state solution of the 
corresponding system showing the asymptotic stability of the 
equilibrium state under the condition (3.48), The condition 
(3.48) is similar to (3.28) and in the special case of equal 
convection and dispersion these are same. 


Further, even if (3.48) is not satisfied the equilibrium 
state may still be stable under (3.28). 

(ii) Semi infinite nab it at' 


In the case of semi- infinite habitat the solution of 
the system (3.41) is found by using Laplace transform with 
respect to t and Fourier sine transform with respect to x and 
y and then using the transformation (3.40), the final solution 
of the system (3.4) can be obtained as follows : 
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U. 


"l (x,y ,p ) = VZ-KZ'T e3C J?<sr r x) J 


'ly 2y 


lx m=l *- 


A 


1m 


(p) exp (-a ^ 2 y) 


+ B lm (p) exp (-q^y) 


Vm 


U, 


sin (SS) 


oo , 


(3.50) 

1/2 


w 2 (x,y,p) = ex P x) j_A Sm (p) exp (-g.^ y) 


+ B 2m (p) exp (-ctg/y) sin (3j2) 


where 


A mh>> = 


2 2 
m n 


K (p)6p( m P : i_y(qi m ;D 2y“P- il -- ;D 2 X 


) + f D 2y P(p) e g (m) 


^Im “ q -2m 


2 2 
m % 




N(p)e 1 (m)D ly (p+bfDg x 2-2- _ q 2m D 2y ) - f I> 2y P(p) 0 2 (m) 


Q-lm " ^2m 


m 2 2 

m % 


(3.51-a) 


= 


S 2m<P> = 


2(p)0g(» 1 )l | 2y(<ll m D 3y-P -e_:D lx - l 2 } + 8 D ly e l^ m - ) 

Q-lm “ q 2m 
2 2 

PCpJegCmjDgyCp+e+D^ 2_|_ _ - g H(p) e^m) 


where 


9-lm “ ^2m 


(3.51-b) 


™ 2 2 
m % 


Z±m’<l2m = 2S^%~ \} p+e+J) ±x ^f’ )D 2y +( P +lH ‘ I) 2x L 2 } D ly 

± |(p+e+D lx ^f-)3) 2y +(p+h+D 2x ^§~) B ly > 2 (3.52) 

2 2 2 2 -i l/2-i 

-4 {(p+e+D^ ~|-)(p+M 2x -*-) D l7 D 2y + fg } J J 
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e-j/m) 

e 2 (m) 



exp 


exp 


C- &r~ *> Bta <Tr) to 

-IX 

/ ^2 \ . -■ mux \ , 

(- 2 i) — x) sin dx 

2x 


(3.53) 


and N(p), P(p) denote the Laplace transforms N(t) and P(t) 
respectively. 


To investigate the asymptotic behaviour of v^(x,y,t)£ 
i = 1,2 as t - », we use the following property of Laplace 
transform £ See Le Page (1961) p. 315 3 


Lt v-(x,y,t) = Lt p v-(x,y,p) 
t - 00 p - 0 


(3.54) 


Using (3.54) in (3.50) we get 

2 r u l 

v 1 (x,y,t) = l exp 


V. 


X + 


'lx 2D ly 


] 


00 

Z 

m=l 


[ 


-^) + f D zy P Q 8 g W 


exp 




o c 

q lm ~ ^2m p g 

lu ” 0 


ViC“>V m 2x Tf - WV'^VoV^ 


exp 




o o 
3-lrn “ ^2m 


. ,mitxx 
S3J1 


(3.55-a) 



6 5 


v s (x,y,t) = exp [§ 


2D 8x X + 4>y a 

2 2 
m 7c 


- t£ 0 + 6 M c 6 l< m > 

Z - — — 

m=l L 


o n 
Q-lm ~ ^-2m 


oxp [i? m )1/2 y_ 


P o 0 2^ m ^ I) 2y^ e+I) lx ”" T 2 " ^2m ^ly ^ “ s I) ly N o e i^ m ^ 


o o 
^lm ~ q 2m 


exp 


y] sin ( 


nmx\ 

T“V 


(3 ,55—b) 


which would exist provided 


IT 2 V 2 Tl 2 tt 2 tt 2 -n 2 

U V D ta n U 2 V 8 D^n .2 , _ 

(zfi — + ITT— + — o”“'vZT5 — + 2T) — + t — J - s > 0 ( o . 56 y 

4xJ lx 4JJ ly L 2 J - U 2x 4IJ 2y l 2 


where q^ m ? i = 1,2 is calculated from (3.52) at p = 0. 

For eq.ua! convection and dispersion the condition for 
stability (3,56) is the same as (3.37). Further, even if the 
condition (3,56) is not satisfied the equilibrium state may be 
stable under condition (3.29) also. 


3.4 COMPETITION MODEL UNDER ZERO LONGITUDINAL DISPERSION 

In the previous article, we have studied the linear 
stability of competition model for nonzero dispersion coefficients. 
However, it may be noted here that the transformation (3.14) and 
the solutions thus obtained are not valid for D 1 (u 1 ,u g ) = 0 or 
Dg(Ui,u g ) = 0 when ^ / 0 or U g / 0. Therefore, to, understand 
the behaviour of the solution insuch cases, in the following we 
discuss the stability of the system when D^ = 0 or Dg = 


0 



This case is applicable when the dispersive migration in a 
particular direction is zero due to dense population of the 
species arising out of non-homogeneity of habitat or wnen the 
convective migration is dominant over dispersion in that 
particular direction. 


In the following the stability of equilibrium state is 
discussed for the case of eq.ual convection and dispersion for 
two species when D lx = D gx = 0, 

(i) Finite habitat 


Using the transformation 


X = exp 
X^ = exp 


[ 


vy v 2 t 
2D " 4ir 

v 2 t V y 
4D “ 2D 


] 

] 



(3.57) 


tbe solution of the system (3,4) is obtained under initial and 
boundary conditions (3.5) - (3,10) as 


oo 

v i (x,y,t) = exp(§j) 2 ^ m G^(t)sin (SSL); t < x/U (3.58) 

OO 

v i (x,y,t) = — exp (|J; 2^ m S£(t) sin (5jE£);t > x/U (3.59) 


for i = 1,2 
where 



67 





a 

1/u/ 

o 


V i;L (T/U) exp 


Dm 2 ir 2 V 2 
" l2 4D * 



(3.61) 

(3.62) 


We discuss the following two cases i 

Case Is U = 0# In this case the population densities 
v i (x,y,t), i = 1,2 is given by (3.58) and, as before, its limit 
as t — °° can exist provided 

P^(m) - s > 0 (3.63 ) 

is satisfied and can be obtained as 


v“(x,y) = exp (!$) 


rVy 


T|(x,y) = ^ ex P (g) 


2^ m F-J^m) j^Pj_(m) H 0 -fE 0 ^| sin (SgL) 


(3.64) 


OO mni 

2 nfl^m) [^P^(m)P 0 -fH 0 J sin (SgE) 


m=l 


where 


■^i 1 ( m ) = - s 2 ] 


It can be noted that (3.64) is also the steady state solution 
of the corresponding system showing the asymptotic stability 
of the equilibrium state under the condition 


V s . Hit® 


(3.65) 


s > 0 
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Comparing the conditions (5.28) and (3.65) it is also seen that 
the condition (3,65) is a one dimensional version of (3.28), 

Thus, it can be concluded that the epu'' 1 "brium state which is 
stable in two dimensional habitat may not be stable in one 
dimensional habitat, 

Further, even if the condition (3,65) is not satisfied 
the equilibrium state is stable under the condition (3,29) and 
the limiting solution obtained as t - °° is the same as the steady 
state solution of the system, given by 

oo 

v®(x,y) = exp (Jyjj y) Z m Fj^m) sin (2ffi) 

Jj m=l 

(3.66) 

v|(x,y) = exp (|jy y) ^ m £ 0 ^iO) sin (^-) 
where 

F| 1 (m) = |”p^(m) + s] . 

If none of the conditions (3.65) and (3.29) is satisfied then 
the equilibrium state is unstable. 

Case 2 : U ^ 0. In this case the distribution of species as 
t - oo is given by (3.59) and it is the same as the steady state 
solution of the corresponding system showing that, in the absence 
of longitudinal dispersion, the equilibrium state is always 
stable with convection in that direction. 
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(ii) Semi-infinite habitat 


In the case of infinite habitat the distribution of the 
species are obtained as follows : 


v (x,y,t) = / 7^(1) $(x,y,I) ai, t < x/U 

0 

X 

Y i (x,y,t) = / V il (Z/U 1 ) (x,y ,Z) dZ, t > x/U 


(3.67) 

(3.68) 


where 


♦ (x,y,T) = 1“ 3 / 2 

2f nJT 


i = 1,2 


r v 


exp [__2U y " 4^ T 


f(x,y,Z) = 


y 


-3/2 


„ , r «U " exp j 2D y “ 4DU “ 4DZ 

^ V xj 


[ 


L-v- O - O' 


(3.69) 


(3.70) 


In the case U = 0, the finite limits of v^Vg, given by (3.67), 
as t -*■ oo can be obtained as 

.2 


v i(y) = H 0 [ -( X2 + n )1/2 


(3.71) 


T f(y> = p 0 exp [- <^g + ! )1/z y ] 


provided condition (3.29) is satisfied. Further, the distribution 
(3.71) is the same as the steady state solution of the system 
(3.4) under (3.29) showing the asymptotic stability of the 
equilibrium state. It may be noted that the condition (3.29) 
is also needed for the existence of the steady state solution 
of the system (3.4) under boundary conditions (3.6) - (3.10) ► 
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Comparing this case (D^ = 0), with 0 in infinite 

habitat case, it can be seen that (3. 29) is not needed for the 
existence of the steady state solution, Further, the case 
(i.e, D^ = 0), stability of the equilibrium state (in the 
absence of longitudinal convection) remains unaffected by the 
presence of transverse convection and dispersion while for 

^ 0, the transverse convection and dispersion have stabilizing 
effect. 

Again, for f 0, it can be noted from (3,68) that the 
equilibrium state is always asymptotically stable to the steady 
state solution, 

3,5 P£EI PREDATOR MODEL 

In the case of prey predator model, the interaction 
function f^, i = 1,2 is prescribed as follows : 

f l = u l^ a l “ c l u 2 ^ 

(3.72) 

f 2 = ^ g ^ 

The linearized system for the evolution of two species would be 
given by (3.4) with 

b u = 0 

b 12 = “ f < 0 (3.73) 

b 21 = S > 0 

The solution of the system (3,4) with (3.73) under 
conditions (3.5) - (3.10) in both finite and infinite habitats 
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•for = 0 or 1)^ 0, are given by respective solutions in the 

2 2 

case of competition model where s is replaced by -s . 

Noting (3i73) it can be s^-on from the solutions of 
the system in various cases that the equilibrium state which 
is otherwise neutrally stable is always asymptotically stable 
with convection and dispersion giving ncnuniform spatial patterns. 

3*6 CONCLUSIONS 

The linear stability of two interacting species system 
in two dimensional habitats has been investigated by considering 
dispersive and convective migration of the species under time 
dependent boundary conditions. It has been shown that condition 
for stability of the equilibrium state under time dependent 
boundary conditions is the same as for constant boundary 
conditions if (3.9) is satisfied. 

The following conclusions may be drawn from the analysis 
for the competition model. 

1. In the case of a infinite habitat for D^» Dg nonzero, it 
has been shown that the equilibrium state may become stable 
with dispersion and convection in a two dimensional habitat 
provided (3.28) is satisfied. 

2. Even if the condition (3.28) is not satisfied the 
equilibrium state is stable under the condition (3.29) 
involving the boundary conditions. 
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3. The equilibrium state which is stable under the condition 
(3.28) in the case cf finite habitat may not be stable in the 
case of semi- infinite habitat because (3.37) may not be 
satisfied. 

4. The equilibrium state is always stable with ^ 0 and 
= 0 in both finite and semi- inf inite habitats. 

5. It is noted from condition (3.28) and (3.65) that the 
equilibrium state which is stable in one dimensional finite 
habitat would stablize further in two dimensional finite habitat 
due to convection and dispersion in both longitudinal and 
transverse directions. 

6. In semi- inf inite habitat for = 0 it has been 

shown that stability of the equilibrium state is unaffected by 
transverse convection and dispersion. 

In the case of prey predator model the equilibrium 
state is always stable with dispersion and convection in all 
cases. 
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CHAPTER IV 


EEEECTS OP DISPERSION ON THE LINEAR STABILITY OP TWO 
SPECIES SYSTEM H\i A T vifO-DIMEN S I ORAL PATCHY HABITAT 

4.1 INTRODUCTION 

Though the mathematical modelling of interacting populations 
in a habitat has become quite sound after the work of Lotka and 
Volterra, it is only recently that the environmental and ecological 
effects have been considered on their evolutionary processes 
( Skellam, 1951 ; Segal and Jackson, 1972} Hadelar , 1974 $ 
Gopalsamy , 1977 5 Kerner , 1959 ) . Generally these effects are 

taken into account by identifying species migration with dispersive 
and convective processes which might have continuously varying 
properties due to environmental and ecological gradients in the 
habitat (Levin, 1974 $ Comins and Blatt, 1971 ). Although 

it has been pointed out by many workers (Huf faker, 1958, 1963; 
Comins and Blatt, 1974*, Levin, 1974 ° 9 Kawasaki et al, 1979) 
that migration of species may depend upon densities as well as 
location of the species, little attempt has been made to study 
such effects on the stability of population models. In such a 
case, not only the interaction terms but also the migration terms 
become nonlinear in the model and it is difficult to study such 
a system* 

One approach to study such a problem is to divide the 
habitat into finitely many sub-habitats (patches) where the 
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dispersion coefficients may be assumed to be constants with 
values different in different patches. Shis method, though 
not exact, can give fairly close approximations e^s regards to 
the effects of dispersive migration on the evolution and 
coexistence of the species (Levin, 1974? Kawasaki et al, 1979). 

In this chapter therefore, we study the effects of 
dispersive migration on the linear stability of two interacting 
species system in a two dimensional patchy habitat by following 
the above mentioned approach. Though the space dependence of 
interaction rates in the model would give greater realism, it 
would undully complicate the interpretation of the results and 
hence in this and in the following chapters these parameters 
have been assumed to be the same in ail patches (Comins and 
Blatt, 1974). To see the importance of dispersive migration, 
the particular case of Vol terra type competition model has been 
investigated in detail as the equilibrium state in this case is 
inherently unstable . 

4 .2 BASIC EQUATIONS 

Consider the evolution of two interacting and migrating 
species in a two-dimensional heterogeneous habitat [0 < x < L, 

0 < y < Bj • Identifying the migration of the species, by 
dispersion process depending upon the density gradients of the 
species, the governing system for the evolution of species can 
be written as. 
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du. 


a au. a 

f i (u l' u 8 ) + H LD lx (x,y) + fy & iy (x,y) 


au. 
ay - 1 


(4.1) 


where u^(x,y,t) is the density of the ith species at any 
point (x,y) and at any instant t. The terms f^(u 1 ,Ug), i=l>2 
are due to the interaction of the species end (D. ,D. ) denotes 
the dispersion coefficient of the ith species and may vary with 
sp ace . 


As pointed out earlier, divide the habitat into finitely 
many patches, say !, such that the dispersion in each patch is 

constant, let jth patch be defined as 1. 1 < x < 1., 3 = 1 , 2 ,...! 

3““ 3 

with L = 0 and 1^= 1. If u. ., (D- .. ,D. . ) denote the 
density and dispersion coefficient of the ith species in the jth 
patch, the evolution of ith species in jth patch can then be 
written, from (4.1), as 


Qu. . 

rr-r^- = f . 


at 


1 D 


+ D- 4 


xox 


a s u. . 

D 

ax^ 


fin 

W ay 2 


L . 1 < x < L . 

2- 1 0 


1=1,2,; 3=1,2,...! • (4.2) 


Since it has been assumed that the interaction of the 
species remains unaffected by the patchiness of the habitat 
the nontrivial positive equilibrium state (u^.., u^.. ) of 
the above system can be obtained by solving 

f = 0 ; i=l,2, ; 3=1,2,...! 

1 J 

which are the same in all the patches i.e. u^- = in . 
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The system (4.2) can be linearised by writing 

u. . = u. + v. . ; i= 1 , 2 5 3=1,2,,..,! 
x ^ i 13 

and neglecting square and higher order terras to get 


dv 2 a 2 v., 6 2 v, , 

~c— "■ = Z b . . v, .+ D . . — - — ■■4 + D. . — 

dt t Ss1 ik k 3 i 3 x Sv 2 107 


dx ( 


ay 


*3 ) 


where 

sf ii 

^ik = 

0 U k3 Ul ,u 2 

and v • . , i=l,2 5 3=1,2,...! are the small perturbations about 

J- J 

the equilibrium state. 

The following initial, boundary and matcning conditions 
may be associated with the above system. 

Initial conditions : 


v i; j(x,y,0) = G^(x) 51, < x < L.. (4.4-a) 

4 • 

such that 

0 «i-i)<V = W 

(4.4-b) 

The equation (4.4-b) gives the matching of initial conditions 
at the jth interface. 

Boundary conditions : The boundary conditions may be any of 
the following types: 
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1. Reservoir conditions 


v il (0,y,t) = 0 
v iK (L,y,t) = o 
^(x,0,t) = 0 
v i ^(x,B,t) = 0 

2. Flux conditions: 


i=l,2 
3= 1 , 2 , t * 3 


( 4 . 5-a ) 


3 . 


or 


6Y il 

ax 


= o 


av. 


i! 


dx 


= 0 


at x = 0 

at x = L « i= 1 , 2 , $ 3=1,2..,! (4.5-b) 


__il = 0 at y = 0, 3 

Mixed conditions : 

v i 1 (0,y ,t) = 0 

>y » "k ) = 0 

8v ii 

= ° y =0 > B for 3=1,2....! 


( 4 • 5— c ) 


dv 


il 


ax 


0 at x=0 


ilil 

ax 


0 at x=l 


v i; ^(x,0,t) = 0 
Vi^(x,B,t) = 0 


i=l ,2 


3—1 , 2 , . . »! 


(4.5-d) 
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Matching conditions s Keeping in view of the continuity of 

-j- 4- 

species densities and fluxes at the interface of (j-l) and 3 
patches 3=2, 3,.. .N the matching conditions are 


av iM-l) av ii 

D. / . 1 x = d. . — 12 a t %=& . 

1(3-! )x Qx 13X ax a 3 


(4.6) 

(4.7) 


for i=l , 2 , 5 3=2, 3,... N. 


how in the following the linear stability of Yolterra 
type competition and prey predator interaction models have 
been discussed. 


4.3 COMPETITION MODEL 


The volterra interaction function for the competing 
species is given as follows: 

f ^ ( U^ j Ug ) ^•■j^( Uj^”*C .^Ug ) 

f 2^ U l ,U ’g) “ U 2^ a l“ C 2 U 'l j) 

and the linearized system corresponding to (4.2) can be 
obtained from (4.3) as follows: 


(4.8) 


a 2 v, . 


a 2 

a v 


rr" 1 - - fv 23 + D ljx + D iiy ^2^ ; J '3-l <x - 1 j’ 


3t„ . 


aT- = -S T 1J + D 2jx -r' 1 + D 2jy 


3 %i 


!!zn 


3— 1,2,»»»N (4.9) 


dx~ 3y‘ 


where f = c^ag/Cg 5 g = c 2 a i/ c i' 



It may be noted that these equations are of reaction 
diffusion type which are linear and coupled, io reduce them 
into diffusion equation, the following transformation is used. 


T ij 


_ V 2 i 


-sinh (st) -1/r cosh(st) 


1 

i 

r^l 

1 1 
J 1 

? & 
ro 

CJ. 

I 


j f , 2 , . . 

(4.10) 


where 


r = Yf?g i s = fig 

The system (4.9), in transformed form, can be written as, 


_ a 2 5 2 

8 z ±i a z ± . a z ± . 




ay' 


L . .,<x < L.. 

1 — J 


(4.11) 


The transformed initial and matching conditions are 
z 1; j(x,y,0) = rG 1 (x) 


z 2j( x,y ,0 ) " - G p( x V r 


(4.12) 


and 



(4.13) 

dz. . 

= V. . r—— 1 “■ at x=L . 

13X 8x 3 

(4.14) 


dz . / ■ - 
r> r(.i-- 
D i(j-i)x ax 

while the transformed boundary conditions would be given by 


(4.5) with v . - replaced by z± i 

X J X J 
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4.3.1 DISCUSSION OB STABILITY UNDER RESERVOIR CONDITIONS 


Solving the system (4.11) under zero reservoir Boundary 
conditions, initial conditions (4.12) and matching conditions 
(4,13) and (4.14) we get 


oo 00 


z 13 (x,y,t) = c m n p mn^ x} sin ^ B } exp (“\nn t} 


( 4* 15 ) 


oo oo ^ 

z 23 (x,y,t) = s £ ^nW 30 sin (2 S X) ex P ( - A A t} (4<16) 

^ m=l n=l ^ 


where 


I 'mnl (x) = 


sin( ^) 
fDllx 

sin( isSS±l ) 


^ a llx 

^■mni — 1 ^ 


(4*17— a) 


W X) = C0S ^tnn(i-l)( L o-l } ( 4 - 17 " p ) 

« SJ '\ i- 


fHmik- sln r^ S^Tl j=Rl 


+ iPW l 


ljx A mnj 


. / ^mnl x 't 
sin( — ) 


W 1 * = 


Yb 


2 lx 


si n(i2Sid) 

YD gl x 


Y-Dl 


d 






(4,18-a) 



81 


W x) “ 003 (l^) 

2 jx 


+ Sin [-iffijifihli. 


^ D 2 jx x mn;j 


TTD 0 , 


2jx 


dx t I 'in( j_i)( x )D 


x=L 


D- 1 


nmj 


hnnj 


X mn ” D ijy 

- D 


2 2 
n it 

B 2 

„2 2 

n it 


ion 2 jy b2 


L B 


o _ 0 0 

°mn 


/ / z 1 (x,y, 0 ) ! m (j) sln( 2 |i) dxdy 


ll I ' mn (x) sin <tF>II 


G ’ - 0 0 

mn 


I f z g (x,y, 0 ) ^(x) sin ( 2 g£) dxdy 

n n ' 


1 1 0=0 sln(£g£)|| 


(4.18-b) 


(4.19) 


(4.20-a) 


(4.20-b) 


L B 

| | A || 2 = / / A 2 dx dy 
o o 

and ^mn* X mn are r00 " fcs °f ’ tile following equations, obtained 
by using zero boundary conditions at the end x=L in (4.16) with 
(4.18-b) for 3 = 1 . 


J l(M-l)x 

D 1Bx 


x mnl ^mnfN-l)^-!^ 

S' 


cot 


rani 


(1—1 


'H-l^ 




3=0 


(4.21— b) 
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D 


2(U-l)x 'SanU 


D 


2Nx 


^ ^2Nx dx S'mn(N-l)^ x ^3 x= L s - 1 


. cot 


[ 


X mr\w( L ~ L i\r-i ^ 


= c 


^ D 2Bx 

(4.21-b) 

respectively. Substitution of (4.21-a) in (4.17-b) for j = N 
gives, 

g mn(U-l)(%.l) r X mnN^ x " L N-l^' 


^mn]tf( x ) 


sin 




x mnE^ L “%-i ^ 


— 


. sin 


V* 


lUx 


(4.22) 


V^asx 

Similar expression for F ^^ x ) is also obtained. 

Using the transformation (4.10), the final solution of 
the original system (4.9) can be obtained as follows : 


°° oo p 0 

lj (x,y ’ t) = Jji L°™ *W x) r 00sb (st} exp 


V 


+ °iiii 3 W x) slnl1 (st)exp(-X^ sln 


(4.23) 


00 oo _ 2 

v 23 (x,y,t) = 2 2 ¥ xnn^ } sinh (si:} exp C-*nn t} 


m=l n=l 


+ C >;»,(*> «p ^] 8 “ ^ 


It can be noted from (4.23) that the finite limit of 
v..; i = 1,2; j = 1,2,...,N as t - would exist provided 

4a" 8>0 

(4.24) 


\ x l - s > 0 
mn 
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which are combined to give the following condition for stability 


X ll X 11 “ s 2 > 0 (4.25) 

where denote the minimum values of 

respectively. It can be noted that zero values of Aj^ 

are also the solution of (4*21) giving trivial solution in 
(4.23) and hence xj^ are taken to be nonzero* Thus, 

the equilibrium state under the condition (4,25) is asymptotically 
stable to the steady state solution (zoro solution). 

It can also be noted from (4.23) that the condition 
(4.25) is common for all patches and thus stability in a 
particular patch not only depends upon the dispersion coefficient 
and length of this patch but on other patches also. 

To see the effects of dispersion on the stability of 
the equilibrium state more precisely, we consider the particular 
cases of habitats with two and three patches only. In case 
of two patches the equations (4*21) determining A^ and . 

reduce to the following form 


mill 


oot + y 

fc- /-Her-.*— — I— 4 » 


/S 


mnl 


cot E 


■Y*nx 

y A*'**"’ 


12x , 

T3 x 

J llx 




2 + 


2 lx 


oot C ■ — 3=0 (4.26-a) 

W>12 x 

5 -is oot (4.26-bj 

D 21x ml2 
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■ B ' or e l u al disper c,i 

it o-vn -hp. ^ 011 00effici ents in y directions, 

it can be noted from (a 1q x x 

that 


mnl 


W)= X 


2 

mi 

• 2 


X c - T 2 
mn 2 ~ Vj 


,= r 


,2 


«n2 mn 
reducing the e, uatlon ^ ^ 


cot r ^n L i - /5~r x fL-L x 

L Y^- 1 ! + T btt «<* C-^Lr- >° (4.2V) 

11 x Y^ 

and similar equation fo „ 

. . %r, can also be written from 

(4.26-b). ^ 


In case of three 


dispersion co e f fi C i ents 


patches, the equation (4.21) for equal 


111 y directions, gives 


'E 


V Tr^ 2 - + o.n+ r- X mr.( L “ L o) 


^ + 00t D 


YB 


1 • 


13x 


T^llx «»»E ^=kl sinC slr Imnh 

^jl2x 


p ~T* rr — ^—r 

c °sL,Tizr:Jco s r WzcS -i - 

llx ! “ 




-HMWM 


^Ux ** VT3 


— r . • i— Xmn^l — i . i— Xnin^9“^1 ^ 

J -f D i 2x sin C >1* C - - 


I2x 


fD 


Hx 




= 0 


and similar equation f nv . - , 

Ui \jjjj can also obtained from 


(4.28) 


Equation (4.27) is so1vq , 

v ea numerically for the minimum root 

and it •? « t 

plotted in figure 1 for various values 


x u of X 


urn 
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rs 


— r 1 

of Dg = “7*" &nd D j_i x » L-jA "being fixed. It is clear 

from figure 1 that X^ iiicreases with Dg for all values of 
1^/1 and this increase in X^ is enhanced as L-j/t decreases 
for Dg > 1 "but reverse is the case when Dg < 1 , 

Again equation (4 .£8) is solved numerically said the 

r flttZ 

graph is plotted in figure 2 between A. 1 and D, = V 

11 6 

from which it can he noticed that X^-. increases with 
D^ f°r fixed Dg, L^/L, lg/L. Further, the increase in A^l 
enhances with Lg/L for D^/Dg >1 for fixed L^/L while 
opposite is the case for D^/Dg < 1 • 

By induction, it can be argued that these results are 
valid for habitats with finitely many patches also. 


Thus, keeping in view of the condition (4.25) and 
expression (4.23) it can be seen that convergence of unsteady 
state solution to the zero steady state solution increases as 
the dispersion in the patchy habitat increases. In other 
words the degree of the stability of the equilibrium state 
increases with dispersion. 


4.3.2 DISCUSSION OF STABILITY UNDER FLUX CONDITIONS 


As before in the case of flux conditions, the solution 
of the system (4,9) can be obtained as follows 
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v ld U,y,t) 


m=0 n=0 *3 (X) r C0s ^ st) ex ^~ X L^ 


+ ^mW 30 Sinh(st) ex P(-O):° 0s (^) 


» 2_ 
mn 


B 


CO 03 


v 2j (x,y,t) = y e 2 p(-C,t) 


m=0 n=0 ^ 1211 


Torn 


+ C mn F mni( x ^ — G3 P(“^ t )3 ^(^g 2 -) 


where 


L B 


f f z 1 (x,y,0)cos(^) yxjdxdy 


0 o 


mn 


ll P mn( x) cos (tF)!! 


L B 


/ / z 9 (x,y ,0) cos(^) By(x)dxdy 


C ' 
mn 


o o 


IlF^tx) eos (£f£)l| 


COB { 


*mnl x 


W z ^ 


0 


llx 


. .. r X mnl L l, 
cos i ■ • '--•—} 


fE> 


llx 


, Ci x . 

COS { - } 


■^mnl^ 2 ^ = 


0 


2 lx 


, mnl L l , 
cos { - } 


fD 


2 lx 


(4.29) 


(4.30) 


(4.31) 



87 


and aXe given by (4.17-b) and (4.18-b) 

respectively . 

As before, using zero flux boundary conditions at the 
other end x=L, the equations determining x^, can be 

obtained as 


D 


)x _ X mnN ^mn(H~l 

1Nx ^ fe l^mn(N-l)^ x ^ 


:an C 


W L ~W 


^N-l 




1 = 0 


lflx 


P 2(N-l)x _ X nmN ^mnfN-p^N-l^ 

D^iiT A 


2Nx 


^ P 2Ix lx ^mn(N-l/ x ^ 


(4,32-a) 
tan f^g^1 = o 


fs 


2Hx 


x=l 


N-l 


(4.32-b) 


respectively . 


From the solution (4.29) it can again be noted here that 
the stability condition in this case is also given by (4.25) 
where in this case X^, X^ 1 are the minimum roots of (4.31-a) 
and (4*32-b) respectively which could even be zero as they 
would give nontrivial solution in (4*29 )» fhus, when X 11 “ 

^ = 0; the condition (4.25) will never be satisfied and 
hence the equilibrium state which is unstable without dispersion 
remains so with dispersion in patchy habitat under flux boundary 
conditions* 
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4.3.3 DISCUSSIOH OF STABILITY UNDEB iilXED CONDITION'S 

Proceeding in the same manner as before, the solution 
of the system (4.9) under mixed boundary conditions of 
type (4.5c) is given by 

OO GO 

V lD (x, y , t ) = m l ± n l Q &w?wi}W r cosh(st) exp(-A^t) 

+ C mn F mn;j( x ^ sinh(st) exp (- A^t)^} cos (-^) 

(4.33) 

OO OO 

v 2d (x ’ y,t) = slnh(st) ^(-4^ 

+ °iiiii i 'imT X ^ C ° S r <St ^ ex P(" x 'jtn t ^ cos (^) 

where C^, are given by (4.20) and F^, F^ are 

given by (4.17), (4.18) respectively. Also A , in 

this case are the roots of equations (4.21). From the 
solution (4.33) it is obvious that the stability condition 
in this case also is given by (4.25) where An, A^ , the 
minimum values of A^, correspond to nt=l, n=0, 

[see(4.19)3 . As discussed before A , Aj^ remain nonzero 
in this case and would not depend upon the component of 
dispersion coefficients of the two species in y direction. 

As these minimum roots increases with dispersion of the 
species in x direction enhancing the degree of stability of 
the equilibrium state# 
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However, in the case of mixed boundary conditions of the 
type (4.5 d), the solution can be obtained as 


oo oo 

v lD (X,y,t) = iJo n=l ^ 0mnFim D (x) r COsh(st) ex P(“Vn t) 
+ G mn F mn^ X ' 1 sinh (st) ex P(“ sin 

oo oc 



(4.34) 


Where C^, 0^ are the same as (4.30), S^, are given 

by (4.31) and $ m y JjJmj *7 (4.17b), (4.18-b) respectively. 
In this case X^, ^ are the roots of equation (4.32) 

which could be zero also. Thus the minimum values of 


X* 


mn 

and would not depend upon the dispersion coefficient in x 
direction [see ( 4 - 19 ) H • As the c^i^io 11 for stability in 
this case is also given by (4,25) it can be noted here that 
degree of stability of the equilibrium state increases with 
dispersion coefficient in y direction. 


i.e. x-t-ij x‘ -ii w iH now correspond to nt=0, n=l 

JLJL 


Comparing these two cases of mixed boundary conditions 
with the case of flux boundary conditions, it can be noted 
that dispersion in a two dimensional rectangular habitat with 
reservoir boundary conditions along two opposite sides of the 
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habitat can make the equilibrium state stable which is unstable 
with zero flux boundary conditions all around the boundary. 

4.3.4 THE CASE OE OHS DIMENSIONAL KABIlAl 

In the case of one dimensional habitat the solution of 
the system (4.9) with reservoir conditions can be obtained 
as follows : 


Y lj “ C B n r =°sh( s t) exp (-w B t) 

+ B n T nj^ x ^ sin3l (st) exp (-w^ 2 t)] 

OO 

- E -£B n 3 n j ( x ) sinh(st) exp (-w 2 t) (4.35) 

+ T^(x)/r cosh(st) exp (-w^ 2 t)] 


where 


3 = 1,2, ...K 


sin ( 


w x 

S— ) 


= 


^llx 


W L 1 

sin ( -- ) 


(4 ,36-a) 


fD lX 3 


w ( x-1 . 1 ) 

T nj (x) = G0S C ' x 'n( j-1 ) (L d-1 ) 




(j-D V A 


(4.36-b) 


GL( 3 -i)W1_ : 


2jx 


n 


^lox 


X= V 
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w! x 

Sin n 


^l(x) = 


YD 


2 lx 


w* L-, 

sin C™=zH 


fD, 


2 lx 


(4.37-a) 


^nj^ x ^ cos C 


H A! x - L .i-:9 


^Sjx 


■1 


(4.37-b) 


+ ^SU=iii_ sin a 


^Sjx w i 


f B 9-i 


1 dx Bn(d-1)^-I 


2jx 


^ L j-i 


B = 
n 


/ Z, (x,0) 1 (x) dx 

o • L 


(4.38-a) 


2, 


n 


B 


/ Z i (x,0) T^(x) dx 

f _ o __ 


n 


rn I I I 2 


n 


(4.38-b) 


w • W* can be determined from the following equation 
n* n 


D l«l-l)x w n T n(M-l) £a=lL 


D 


cot {3 


lNx 


S l’r-(iI-l)^ X Y-= 




3= o 


X=L, 


H - 1 (4.39) 

D gfM-l)x W n <W cot r W ° (1 ~^ -l = 0 


D 


2Hx 


f%x dx 


^ I) 2Kx 


K-l 
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She condition for stability of the equilibrium state 
in the patchy habitat is obtained as 

W n W A 2 *“ s 2 > 0 (4.40) 

As seen before, the minimum roots w^, increase 
with dispersion showing the enhancement of the degree of 
stability of the equilibrium state in this case also# 

Further, comparing X^, xj^ with w f in view of 
(4,19) it can be noted that even if the condition (4,25) is 
satisfied, the condition (4.40) may not be satisfied. Hence 
an equilibrium state can be stable in a two dimensional patchy 
habitat even if it is unstable in one dimensional patchy 
habitat . 

In case of flux boundary conditions the solution of the 
system (4.9) is obtained as (4.35) with 
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and 3=2, . . .35T axe the same as (4.36-b), (4.37-b) 

respectively • The equations determining w n , can he 
obtained as 


• D 1 (E-l)x 
5 1 Nx 


D 2(N-l)3C 

^gNx 


W n T n(E-l) (%-j) 
^" S llx lx &n(N-l)^ x ^ 


w (L-l x ) 

, tan f — — ^ - 3=0 


V* 


INx 


x=L 


3S— 1 


(4.42) 


w n T n(N-l)( L E-l^ 

^2Nx lx Ki^-l)^^ 


W (L-Iur ) 

tail r - n - " - -g = ^ -l= o 


f D 


2Nx 


x=L. 


N-l 


which give zero minimum roots for the nontrivial solution 
(4.35). As the condition for stability is the same as (4.40) 
in this case also, it can be noted that the equilibrium state 
is unstable under flux boundary conditions. 

Thus, it is concluded that the equilibrium state which is 
always unstable under flux boundary condition in one dimensional 
patchy habitat may become stable under mixed boundary conditions 
of type (4 .5-d ) in a two dimensional habitat. 


4.4 PREY" PREDATOR MODEL 

In case of prey-predator model the interaction functions 
are. prescribed as follows 

£ ij = u ij (a i - “i'V 
f Sj = U 2J*“V 0 2 U 13 ) 


(4.43) 
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The Uneari zed 


ca21 be obta 


s yst< 


3-in.ed 


’ ln case, corresponding to (4.2) 


- 0a (4.8) 


^ v ^ . 

~d^ = -a 


T&- 


83 + »■, . Idii 8 v li 

J l3x + D„ . hi 

Qx 2 Ijy av 2 




13 + . ^2.1 


Using the foliowj. 


P + T) u 2.i 
dx 2 2 3y 3y2 


^ tr ^sforniation 
COs ( s *) sin (at) 

“** <st) -1/r cos(st) 


system (4 , \ 

' *^) ai 

w e get 


J olvin. 


(4.44) 


(4.45) 


■S under boundary conditions (4.5_ a ) 


®=1 n=i C G ma * F rmi oos(st) exp(-X® t) 


+ 0 * p t o 

^mnj sin ( st) exp(-\^t)J 


v 2i( x »y>t) 


: 1 n=i ^mn^mnj ^ sin(st) exp(~x^t) 


— Q | > 

103:1 r F mnj( x ) cos(st) exp(- x'j^t )]] 

where p Pl sin 

, mn 3’ P mni» 0 B n 

(4.1? )-(4.2o) 103:1 u mn> ane defined, as before, in 

9* 1 ’ *mn can be obtained, from (4.21). 



.Analysing as before, it can be noted from (4.46) that the 
equilibrium state which is stable otherwise remains so even with 
dispersive migration of the species. 

Further, similar analysis can be carried out for the boundary 
conditio 310 (9-b), (9~c), (9-d) from which it can be noted that the 
stability of the equilibrium state increases with dispersion coeffi- 
cients involved with reservoir conditions but would not depend 
up 0 n dispersion coefficients associated with flux boundary 
conditions • 

4 15 CONCLUSIONS 

In this chapter, the linear stability of competition and 
p re y_predator models in two dimensional rectangular patchy habitat - 
has been discussed and the following conclusions are drawn from 
the linear stability analysis. 

q In the case of competition model, the equilibrium state 

which is unstable otherwise may become stable with dispersion 
when reservoir conditions are prescribed on the boundary of the 
habitat and the increase in dispersion coefficients of the 
species enhances the degree of stability of the equilibrium state. 

II The equilibrium state for this competition model remains 
unstable even with dispersion under flux boundary conditions. 

The dispersion in a habitat with reservoir conditions 

X- X 4* # 

along any two opposite sides of the rectangular habitat can 
stabiliz® the equilibrium state for competition model which 



96 


is unstable with zero flux boundary conditions all along the 
b oundary • 

iv. The equilibrium state in a competi ^on model which is 
stable in a two dimensional habitat may not be stable in a 
linear one-dimensional habitat, 

y. For prey predator model, the equilibrium state which 
is otherwise neutrally stable becomes stable for all types 
of boundary conditions discussed here. 
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CHAPTER V 


EPPECTS OP CORVSCfIVE Al'j_D njc-^ ., T - 

Ib PATCHY HAPITaJs-I 1 ^ 1 ^ MIGRATI0K 


5.1 INTRODUCTION 

In Chapter IV, the effects of riio 

1 dls Persiv e migrations on ' 
the linear stability of two interacting 0 

m g species system has been 
investigated in both one and two dimer cn™ i 

nsional finite habitats 

and it has been shown that dispersin-n 

P 31011 caH stabilise an otherwise 

unstable equilibrium state. 

It may be pointed out however, that tv, - + . 

t tnat the migration of 

the species can also he identified by convection ^ 
arise due to biased dispersion (Comins and Blatt, 1974) Ihe 
convective ^ration of the species can ^so arise because of 
sudden changes in ecological and environ ra ental conditions in 

cases of heterogeneous habitats, it ^ v 

can be remarked that 

convective processes might play a dominant mi 

ant role over dispersive 

processes in affecting the evolution sm* . . 

a^u coexistence of the 

species in a given habitat and thus the study of convective 
effects even in absence of dispersion ma y be of interest. 

Keeping the above in view, i n x 

tnis chapter, the linear 

stability of two species system (Voitov.., +. 

^ uirerra type competition and 

prey predator models) in heterogeneous -h^ * - 

6 US two dimensional habitats 
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consisting of finitely many connecting patches (sub-habitats ) 
has been studied by taking into accounts the effects of 
convective migration in both longitudinal and transverse 
directions. The effects of dispersive migration in transverse 
direction has also been investigated. 

Since it has been shown in Chapter II and III that 
the effects of species migration on linear stability of 
equilibrium state in cases of equal and unequal convection 
and dispersion are similar, in this c nap ter we therefore 
consider only the case of equal convection and dispersion 
for the two species wherever it is relevent. 

5.2 SASIC EQUATIONS 

let us consider the evolution of two interacting 
and migrating species in a two dimensional [] 0 < x, 

0 < y < B~[ heterogeneous habitat divided into finitely 
many (say N) patches where jth patch (3 = 1 , 2 ,...,N) is 
defined as £ 1 .; -1 < x < L-, 0<y< B^j with L = 0, Then the 
evolution of two Interacting and migrating species is governed by 
the following system of second order partial differential equations 
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3 u - . Su. 8u. . S 2 u. . 9 2 u. . 

+ u. . + V . . — il s f . . + D - — + D - -ii 

3 t 10 ax 13 3y io + iix + ijy 2 


i = 1 , 2,1 3 = 1 , 2 , . . . ,H. 


(5.1) 


where f. . corresponds to the usual Volterra term for xhe 
X J 

interaction of ith species in jth patch, U. . , V. . are the 

10 10 

convective velocity components and (D. . , D.. ) is the dispersion 

1 JX X Jj 

coefficient for the ith species in jth patch. 


To simplify the analysis, it is assumed that the dispersive 
migration in the x direction is negligible in comparison to the 
dispersive migration in the y direction. This is feasible in 
the case of narrow habitat and when the convective velocities 
are dominant. It is further assumed that the migrative ability 
of the two species are same in respective patches. Keeping these 
in view, we have 


D. . v 0 
ID* ~ 


i = 1,2 


D ioy " D 2jy ~ D o 
u io = u 2o = u o 

V lj = V 2 j = 7 3 


j = 1,2, ... ,N 


(5.2) 


Assuming further that the nature of interaction of the two 
species are same in all the patches then the equilibrium state 
(u 13 = u lf u^ = u 2 ) can be obtained by solving 

f i3 = 0 5 1 = 1,25 3 = * ,N * 
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Linearizing system (5,1) about the equilibrium point (u^, Ug) 
after making use of 


u. . = U- + V. . 
10 1 10 


l = 


■ 5 ^ 


in (5.1) we get 


(5.3) 


9v . . 3v . . dv • . 

+ u, «=3JL + V, 


at 


o 


o ay 


L 

k= 1 


b . , v, . + D . 
ik kj 3 


a 2 v, . 

3a 

o 

ay‘ 


(5.4) 


9f i 

where b ik = (— )_ _ I i = 1,2; j = 1,2, . ..,N 

k u x , u g 

where v. ., i,j = 1,2 is the perturbation density of the ith 

-L J 

species in the jth patch. 

The system (5,4) is associated with the following initial, 
boundary and matching conditions ; 


V i0 

= 0 

? i = 

: 1, 

» 2 5 0 — 

1 ° 1 

-i-jfOy#* » ji 

.■i at 

t = 0 

(5.5) 

y ij 


> 0 

0 

= 1,2,... 

f N at y 

= 0 









(5.6) 

y 2j 


> 0 

0 

= 1,2, .. . 

at y = 

= 0 



•rr> 

l> 

= 0, 

y = 

B 

for i = 1 

, 2 , 0 = 

1,2, 

. .. ,1 

(5.7) 

< 

H* 

H 

= 0, 

at x 

= 

0 for i - 

= 1,2 



(5.8) 

V . . 


v. • 

at 

x = L . 

for i = 

1,2 



1,0 

-1 “ 

i , 0 

o-i 










0 = 

2,3, 

. ,.,N~1. 

(5.9) 


The conditions (5.6) are the reservoir conditions prescribed 
at y = 0 and the conditions (5.9) are the matching conditions 
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prescribed at the interface of jth and (j-l)th patch, j = 2,...h. 


5.3 COMP Eli flOh MODEL 


Considering the case of Volterra type competition model, 


the interaction functions f^, fg^; 3 = 1,2 are given by 


f lj u lj ^ a l “ C 1 J 
f 2j = u 2j (a 2 - c 2 u lj } 


(5.10) 


which will give the positive equilibrium state as (ag/c^a^/c^). 
Using (5.10) the linearized system about the equilibrium state 
can be obtained as (5,4) with 


b u = 0 


b 12 = - f = -°l a 2 / c 2 < 0 

b 21 = "S = _c 2 a l/ c l < 0 

(5.11; 


5.3.1 FIUI'IE HABITAT 


For the finite habitat the solution of the system (5.4) 
with (5.11) in the first patch (0 < x < L^, 0 < y < B) after 
using the conditions (5.5) - (5.9) can be obtained as follows 


T il( x ’ y,t ) 


v il( x ’ y>t ) 


2D-jH; 


2D^tc 

B 2 


£ mG. 1 (t) sin(mxy/B)| fort < x/U. , 

m=l 11 1 

(5.12) 

OO 

£ mG . 1 ( x ) sin(mrcy/B); for t > x/U. , 
m=l 11 1 


(5 . 13 ) 
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where 



a 

_ r 
- J 

0 

1^(1) exp 

rV 

L SD j 

- IV(m) I J dl 





a 



rVw i 



ho (a) 

= i/u^ J 

J 0 

^(l/U^) exp - I'^(m) 1/tT Jdl 

t) 

(5 

.14) 


“ 5 3 

cosh 

(si) - 

P 3 r , 

sinh ( si ) 








(5 

.15) 


= p 3 

cosh 

(si) - 

V 1 

sinh (si) 




B.(m) = 

-p. 2 2 

D .m it 

J 

/B 2 + 

T 2 /C4D 3 ) 

(5 

.16) 


r = 

Yf7i 

» s = 



(5, 

.17) 


Similarly in tne jth patch (h. 1 < x < 1.; 0 < y < B) solution 
of the system (5.4) with conditions (5.5) - (5.9) can be obtained 
as follows 


CO 

v i-?(x,y,t) = 2B,Tt/B 2 E 

J J m=l 

v ij( x >y’ t > = J i 3 + 

^i3 5 

where 


2 

X- 

J,, = 2D.1t/B Z 

13 3 m=1 

mS 13 (- 

f 

-T.y 

J . . = 2/B 2 exp 

13 m=l l 



x-L a 

^•(t) sin (mity/B) for t < — 

3 


X-n 


/ - 


u. 

3 




r. ' 
3 


(x-Lj_!) 


(5.18) 

(5.19) 


(5.20) 


U 3 


^(m) sin (Sjp) 

( 5 . 21 ) 
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K 13 (mj = 


B V-s 

S exp ( — y ) 
o 3 




oosh {S (-U— )} - - u . 

D D 


) r 


x-L 


K 2j (m) 


sinh { s( — -|j“ — ) 
3 

M . 

3 


}^j sin ( nmy/3 ) dy 


(5 . 22-a ) 


B Y -? 


x-L 


3-1. 


id- 


T7 “h 

— j-j -~ 

U 

! 

b 

fs 

9s, 

vH 

1 

*r~D 

u. 

3 

x-L. 

( L 3-l ,y,t 

t 




3-1- 


x-L j_i — | 

inli {s(— jj - — ) } J sin(mmy/B) dy 


(5 .22— b) 


It is noticed from (5.12) and (5.13) that the population 

densities of the two species in the first patch remain unaffected 

by the presence of second patch for all times. However, it is 

clear from (5.18), (5,19) that population distributions in the 

second patch axe dependent on the densities in the first patch 
x-Lf 

for t > — «• — and may increase because of them. It is also 
u 2 

noted here that for, Ug = 0, even the second patch becomes 
isolated for t > 0. 


5.3.2 DISCUSSION OB STABILITY IN TWO PATCHES 

Bor simplicity, we first discuss the stability of the 
equilibrium state in the case when the habitat is consisting 
of only two patches, i.e, N = 2. 
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To see the asymptotic stability and the coexistence of 
two species around tne equilibrium state, we investigate the 
behaviour of v qq(x,y,t) and v ig (x,y,t) for i - 1,2 as t - 

The linear stability of the equilibrium state in the 
first patch is discussed in the following cases : 


Case 1 : Considering the case when there is no longitudinal 
(x-airection) convection in the first patch i.e. = 0, it 
can be noted from (5.12) that as t - v^, v gl tend to 

V li’ V 21 res P ect ively, (Gradshteyn, 1965) where 

Y 


'll' 


p v 1 00 _ f— —i 

(x,y) = 21 )^%/B exp (gg-y) m/D 1 (m) j^I' 1 (m)W 1 - ± sin^^I) 

(5.23) 

Q ~\J~ J OO MUM Mwarf 

v gl (x,y) = 2D 1 Tt/B e xp(g^y) S m/S 1 (m) I , 1 (m)P 1 - gl^Jsin 2|I) 


provided + V^/(4B 1 ) - s > 0 (5.24) 

where D.(m) = P?(m) - s 2 for j = 1 (5.25) 

J J 


It may be noted that (5.23) satisfies the steady state 
form of the system (5.4) with the corresponding boundary 
conditions, showing the asymptotic stability of the equilibrium 
state under the condition (5.24). It can be seen from (5.24) 
that stability increases (convergence is faster) as the 
transverse (y direction) convective velocity increases. Also 
in the absence of transverse convection the stability of 
equilibrium state increases with dispersal. 
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Purtner it can be seen from (5.12) that even if the 
condition (5,24) is not satisfied, the finite limit of v. .,i=l,2 

X -n 

can be obtained under the condition 


Nf = Pfr 


(5.26) 


as follows 

s , , 20 1" Vl 

v 11 (x,y) = 


B 


2 ~ exp (gjp y ) 2 


mN^ 


1 nt=l D-^(m) 


sin (mity/B ) 


(5.2?) 


T f 1 ( x »y) = 


2D^7c v 1 °° mP-j_ 

— r- exp (sg-y) 2 r 

B 2 ^1 m=l D 11 (m) 


. sin (rany /B ) 


where D^m) = m 2 -n; 2 lK/B 2 + V 2 /(4IK) + s, j = 1 


(5.28) 


which is the same as the steady state solution of (5.4) with 
the corresponding boundary conditions showing the asymptotic 
stability of equilibrium state under the condition (5.26). If 
both the conditions are not satisfied the equilibrium state is 
unstable • 

Case 2 i In the case when longitudinal convection is nonzero 
i.e. ^ 0, the density distributions of the two species as 
t -» co a re same as (5.13) which is also the steady state solution 
of the system (5.4) showing that tne equilibrium state is always 
asymptotically stable. It can therefore be concluded here that 
the equilibrium state which is otherwise unstable is always 
stable with longitudinal convection in the absence of dispersion. 
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In the second, patch also we have the following two 
possibilities . 

Case 1 s When Ug = 0, it has already been noted that the 
second patch is isolated and so the conditions for stability 
in this pa.tch are obtained (as discussed in the first patch) 
as follows 

7 t 2 D./B 2 + V 2 /(4D.) - s > 0 (5.29) 

or H. = P.r (5.30) 

J d 

for j=2. 

It can be noted here that the condition for stability (5.24) 
or (5.26) and (5.29) or (5.30) in the two patches for = 0, 
and Ug = 0 respectively are the same as can be obtained for one 
dimensional (y direction) homogeneous habitat. 

Case 2 : Bor Ug ^ 0 the population distribution v^g, i = 1,2 
is given by (5.18) and as t -* °o reduces to 


V i2 = J i2 + J x 


(5.31) 


where Jig, 1 = 1,2 is given in (6.20; which is always finite 


and for = 0 is given by 


J 1 = ID-ic/B" 2 2 exp 

m=l n=l 


r-Bg(m)(—jj— i) ^jn/Dg(n) 


0 (m,n)sin(mity/B) 


J 2 = 


T B 1 (n)V 11 (^i) - fV 21 (^i)"| e(m,n)sin(m7ty/B) 

*— 2 2 -> I 

rz 00 °° r Yp x-L-i -t _ 

^w/B 2 2 exp I g5“ y-Fg(m)(-^ — )Jn/Dg(n) 

111 = 1 11 = 1 

[fl(n)V 2 i(^i)-g V X1 ( e(m,n) sin (muy/B) 


(5.32) 
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provided (5.24) is satisfied and 


oo oo 


=4D 1 ix/B 2 2 exp 

J ' m=l n=l 




8 _ , 




x-L-, 

ll(n)Vii(-^-~) 


J 2 = 4D 1 x/B 2 y-J 2 (n>)(~)]u/B :L i(n)v 2:1 c 


0(m,n) sin (mny/B) (5.33) 

x-L 't —i x— J j • 

-) |n/S 11 (n)7 01 ( 

2 


'1 

TjT^ 


e(m,n) sin (m-jty/B) 

provided (5.26) is valid. The constant e(m,n) is given "by 


6(m,n) = A/2 £{ A 2 + (2=£) 2 * 2 j" 1 {e xp (AB)(-l) m " n - 1 > 

(5.34) 

_ {A 2 + (2 +£ } 2 **}-! {e xp(AB)(-l ) m+n . ^ 

where 

v 1 v 2 

and Y il (T), Ih(n) are defined in (5.15), (5.16) respectively. 

It is noticed that (5.31) is the steady state solution of the 
system (5.4) with corresponding boundary conditions provided 
(5.24) or (5,26) is satisfied indicating the asymptotic stability 
of the equilibrium state under respective conditions. It is 
noted here that the condition for stability of the equilibrium 
state in the second patch is same as in the first patch showing 
that the equilibrium state in the second patch is stable provided 
it is stable in the first patch. However if none of the 
conditions (5.24) and (5.26) is satisfied, the equilibrium state 
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is unstable in both the patches and this situation arises due 
to zero longitudinal convective velocity in the first patch 
and nonzero in the second patch. 

It is seen that when / 0 the expression for 
Jh , i = 1,2 is obtained as follows 

■z 00 00 p X -* L 1 v ? 

J i =4D lV( u i B ) 2 2 exp j -F 2 (m) ( — ) + gg- y 

n K^(n) e(m,n) sin (m-rcy/B) (5.35) 

L 1 x-Li 

where Z i (n) =/ exp { -l’ 1 (n )T } y (-y— + T/U ± ) dT (5.36) 
o 2 

It can again be noted that (5.31) with (,5.35) is the steady state 
solution of the system (5.4) showing that the equilibrium state 
is always stable if the longitudinal convective velocities in 
the two patches remain nonzero. In particular, if the conditions 
throughout the habitat are identical, the equilibrium state is 
always stable in presence of longitudinal convection. 

It may be concluded from the above that the patchiness 
may cause the instability of otherwise stable equilibrium state 
provided longitudinal convective velocity is zero in the first 
patch and nonzero in the second patch. 

5.3.3 SEMI -IBB IN I IE HABITAT 

In the case of semi- inf inite habitat (in y direction) 
the solution of the system (5.4) with (5.10) in the first patch 
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after using the conditions (5.5) - (5.9) can he obtained as 
follows 


t 


v i;] _(x,y,t)=/ 7 il (T) ^(x,y,2)d2 

for t < x/U^ 

(5.37) 

X 

v ±1 (x,y,t) = / V. 1 (~-)h 1 (x,y,z)dz 
o 1 

for t > x/U^ 

(5.38) 

where <(t.(x,y,T) = — — l”^^exp -(4D 

3 2VnD“ L 

J 

• Tf^y-V-T) 2 

(5.39) 

^(x,y,T) Z“ 3/2 exo|- 

3 2f7tD /IT, L 

J t) 

(4D.T/U 3 r 1 (y-V.l/U.) 2 ^] 


(5 = i) 

(5.40) 


and Y^, i = 1,2 is defined in (5.15) for 3=1. 

Similarly in the second patch, the distribution of species 
can be obtained by solving (5.4) with (5.10) and using conditions 
(5.7) - (5.9) as follows 

t x-1 1 

^ig(x,y,t) = / V i2 (l) <frg(x,y,T)dT, for t < -r— (5.41) 

X—*Xj 

v i2 (x,y,t) = E i + E| , for t > (5.42) 

X 

Ef = / ^(Z/Ug) * 2 ( x >y>z) dz (5.43) 

0 

oo 

= f B(y,x) M i (y) dy (5.44) 

o 

where 

t 

M i (y) =/ V i:L (l+5/Y 2 ) * 1 (L 1 ,y,I) dl 


(5.45) 



Ill 


, r 2~ 

V 0 x 


* 1 (L 1 ,y,T) = ^(x^,!) {1-H 

(y-y)X 


/ u 2 

0(y,x) = / 3 exp f - 

if 4xD g x i- 2 U 2 


|" v 2 


E 


exp ( - 


4D 2 x 


___^ } _ exp { - 


(y+y)% 

4DgX 


(5.46 ) 
(5.4?) 

}J (5.48) 


x = x-L 1 , t = t-x/U, 


2* 


and. H is a unit step function. 

It is clear from (5.37), (5.38) that the distribution of 

the species in the first patch are independent of the presence 

of second patch, while distribution in the second patch are 

x-L,, 

dependent on the first patch for t > — . However, for U g = 0, 

2 

even the second patch becomes isolated. 


5.3.4 DISCUSSION OH STABILITY IN TWO PATCHES 


In the first patch, when U^ = 0, it can be seen that the 
finite limit of (5.37) as t - « would exist provided the 
condition (5.26) is satisfied and is given by, 


Y 

v u = N i exp £ y “ ^ 2 + y 

V 21 = P 1 exp y - + y ] 


(5.49) 


It can be verified that the condition (5.26) is also the 
necessary and sufficient condition for the existence of steady 
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state solution given by (5.49) showing the asymptotic stability 
of equilibrium state in the first patch. 

For 0, the solution vn 1 , i = 1,2 tends to (5.38) 

as t -+ 00 which is the same as the steady state solution. Thus, 
it is concluded that with longitudinal convection the equilibrium 
state is always stable. 

In tne second patch for Ug = 0, since this patch is 
isolated, similar results as discussed in the first patch are 
obtained here also and the condition for stability is (5.30). 

When Ug 7^ 0 the limiting population densities as t — <» 
can be obtained as follows : 

v i2 (x,y) = E i + E-j. , i = 1,2 (5.50) 

where E^ is the same as defined earlier in (5,43) and E^ is 
given as follows : 

00 

= f e(y,5) M.(y) dy, i = 1,2 (5.51) 

o 

where 

V — 

M 1 (y) = N-j. exp { y - {p - Qjrip + s) 1 / 2 y/Yf-j_}. 

-L 2 1 

(5 #52) 

9 

V — Y 

5 2 (y) = P x exp { y- + s) 1 / 2 j 

for = 0, provided (5.26) is satisfied and 

Li 

M (y)=/ T.-Cz/U.+x/Ug) * 1 (i,y,z)dz, i = 1,2(5.53) 

x o x 

for U 1 ^ 0 . 
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Thus, it is concluded that for nonzero longitudinal 
convection in the second patch, the equilibrium state is 
always stable provided longitudinal convection is nonzero in 
the first patch. However, if it is zero in the first patch, 
the possibility of instability may arise because (5.26) may 
not be satisfied. 

It can be remarked here that the convection and 
dispersion of the two species has no effect on the linear 
stability of the equilibrium state in case of^ infinite habitat. 
Therefore, in the following, we discuss the general case of N 
patches in finite habitat only. 

5.4 STABILITY IN »N» PATCHES 

The stability of equilibrium state in the jth paten, 
j = 2,...,N is now discussed in the following two cases. 

Case It U. = 0 
J 

Por zero longitudinal convection in jth patch, it can be noted 
from (5.18) that this patch is isolated from all other patches. 
When t -+ oo the limiting populations in this case are given by 
(5.23) or (5.27) for j = 2,....,N under condition (5.29) or 
(5.30) respectively which are also the steady state solution 
of the corresponding system, showing the asymptotic stability 
of equilibrium state in the jth patch. 

Case 2 s U. ^ 0. 

v 

When longitudinal convective velocity is nonzero, it can be 
seen from (5,19) that the densities of the two species in 
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jth patch depends upon the density distribution of the species 
in all the preceeding patches provided IT^. ^ 0 for all k, such 
that 1 < k < j-l. Thus the stability in this patch is dependent 
upon the stability of the equilibrium state in all the preceeding 
patches. 

It can be seer, that by solving the steady state system 
corresponding to (5.4) with (5.10) and using boundary conditions 
(5.5) - (5.9), the solution v®.(x,y), i = 1,2, j = 2,3,.. . ,N 

J 

is given as (5.19) - (5.22) when v- .(x,y,t) is replaced by 

v i; .(x,y) in (5.19), j^i*^ ~ v i;j_i( L ;j_i»y ) 

3 

in (5.22) for j = 2, 3,..., S' and v il (l^__ 1 , y, t) by v^l^^y ) 
in (5.13) for x = 1,2. It may be noted from (5.19) - (5.22) 
that the just mentioned steady state solution i.e. v?.(x,y), 

J- J 

i = 1,2, j = 2,3, ...,H is also the limit of v. .(x,y,t) as 
t - oo, showing that the equilibrium state is always asymptotically 
stable in all the patches j = 2,3,...,Ii with nonzero longitudinal 
convective velocity. In particular, it may be pointed out that 
if the conditions in all the patches become identical the 
equilibrium state would always be stable due to longitudinal 
convective migration of the species. 

Further, it may happen that the longitudinal convective 
velocity is zero in at least one patch preceeding to jth patch, 
say U k = 0, k < j, such that ^ 0 for k < < j. As pointed 

out earlier the kth patch would become isolated from all other 
patches and the density distribution of the species in this patch 
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would be given by (5.18) for j = k and the equilibrium state 
would be stable provided (5.29) or (5.30) is satisfied for 
j = k. However, if none of the conditions (5.29) and (5.30) 
is satisfied the equilibrium state in kth patch could be 
unstable. 

It may be noted that in the jth patch, the distribution 
of species would depend upon the density distribution of 
species in all the patches preceeding to it till the kth patch, 
and would be given oy (5.19) to (5.22) where vu-^Xjy ,t ) is 
given by (5.18) for j = k, i = 1,2. Thus, the equilibrium 
state in jth patch would be stable (unstable) if it is stable 
(unstable) in the kth patch and the condition for stability 
(instability) in the jth patch remains the same as in the kth 
patch. 

From the above discussions it may be noted that 
patchiness in a habitat caused by longitudinal convective 
migration may lead to instability of otherwise stable equilibrium 
state in this habitat. 

5.5 PREY PREDATOR MODEL 

In this case the interaction functions can be written 
as 

f lj = ( a l ~ C 1 ^j ^ u lj 
f 2j = ^"" a 2 + °2 U l3^ U 2j 


(5.54) 
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Using (5.54) the linearized system could be given by (5.4) with 
b ii = 0f i = 1,2 


h l2 = ~ f = " c l a 2/ c 2 < 0 

(5.55) 

D 21 = g = c 2 a l //c l > ° 

Solving (5.4) with (5.54) and using conditions (5.5) - (5.9) 
the solution can be obtained as (5.12) - (5.22) where 


V 13 (I) = N.. cos (si) - P n sin (si) 
Vg^(T) = P^ cos (si) + IL/r sin (si) 


(5.56) 


It can be noted that from the solution with (5.55) that the 
first patch is always isolated but the density distribution 
of species on the jth patch j = 2,3,...,N are very much 
dependent on all the preceeding patches for U^ f- 0, k < j 
but it is independent of all the succeeding patches. It can 
be seen here that in this case the equilibrium state which is 
stable otherwise is always stable with patchiness also. 


5.6 CONCLUSIONS 

In this chapter effects of patchiness on the linear 
stability of two interacting and migrating species have been 
discussed by assuming that the interaction coefficients are 
not affected by the patchiness of the habitat. 
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In the case of competition model in the habitat with 
two patches only it has been shown that for nonzero longitudinal 
convection and zero longitudinal dispersion in both the patches 
the equilibrium state is always stable with or without transverse 
convection ana dispersion in finite or infinite habitats, ihe 
first patch is always isolated but for zero longitudinal 
convection even the second patch becomes isolated and the 
conditions for stability of the equilibrium state have been 
obtained as (5.24) or (5.26) and (5.29) or (5.30) in the two 
patches in the case of finite habitat. It can be noted from 
(5.24) and (5.29) that the stability of equilibrium state 
increases with transverse convection. However, if the transverse 
convection is absent the stability increases due to species 
dispersal. But in case of infinite habitat the transverse 
convection and dispersion have no effect on stability and the 
conditions for stability are (5.26) and (5.30) irrespective 
patches. These results are same as would have been obtained 
in the case of one dimensional habitats. 

For nonzero longitudinal convection in the second patch, 
it is noted that the equilibrium state in this patch is stable 
provided it is stable in the first patch. In particular, when 
longitudinal convection is zero in the first patch, then the 
equilibrium state may become unstable in both the patches, 
indicating that the patchiness may have destabilizing effect 
in both the finite and infinite habitats. 
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For the habitats with 11 patches it has been shown that 
in absence of longitudinal convection in the jth patch the 
condition for stability of the equilibrium state is similar to 
as obtained in the first patch. However if the longitudinal 
convective velocity is nonzero in rhe jth patch then the 
stability (instability) in this patch depends upon the stability 
(instability) of equilibrium state in all previous patches. 

In particular, if longitudinal convective velocity is nonzero 
in all the patches preceeding to jth patch then the equilibrium 
state is always stable in this patch. 

In case of prey predator model, it has been shown that 
the equilibrium state which is stable otherwise, remain stable 
with patchiness also. 
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CHIP TER VI 


EFFECTS OF CONVECTIVE AND DISPERSIVE MIGRATION IN 
PATCH! HABITATS 


6.1 INTRODUC TIOi'j 

In Chapter V, effects of convection on linear stability 
of two interacting and. migrating species system in two 
dimensional habitats consisting of finitely many patches have 
been investigated by considering both transverse and longitudinal 
convection in all the patches and dispersion only in the 
transveise direction (i.e. parallel to the interfaces of the 
patches). However the effect of longitudinal dispersion has 
not been considered. 

Keeping this in view, in this chapter, we study the 
linear stability of two interacting and migrating species system 
by considering various combinations of convection and dispersion 
in patchy habitats. 

6.2 BASIC EQUATIONS 

Consider the evolution of two interacting and migrating 
species in a patchy habitat consisting of two or finitely many 
patches. By talcing into account the convection and dispersion, 
the governing system for the evolution of two species in the 
jth patch can be given by equation (3.4 ) in Chapter III. 

Assuming the interaction coefficients of the species to 
be unaffected by the patchiness of the habitat, the nontrivial 
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uniform positive equilibrium state (u. . = u. , i = 1,2 5 3=1,2,. .. # H) 
of the system ( 3.4 ) can be obtained by solving 



1 3 2 j 2 


1,2 




linearization of the governing system ( 3.4 ) would give 


au ii 
XL + U. 


dt 


ij dx 


dv . . 

__il + v 


8 X 


id ay 


2 

1 b v, . + D. 
k=l lk k 3 


.2 
a v 


iyx 3x 


■ -XL + X) 


ijy 


ay 2 


To study the linear stability of the equilibrium state 
the system (6.1) is solved in one and two dimensional habitats 
under suitable initial and boundary conditions which would be 
specified in particular cases studied. 


It may be noted here that this coupled linear system 
(6.1) cannot be solved in a patchy habitat due to involvement 
of matching conditions at the interfaces of the patches and hence 
the following particular cases, involving some combinations of 
dispersive and convective migration of the species are studied. 

(i) One dimensional patchy habitat consisting of two 

or N-patches with species having dispersion only in the 
last patch while convection in all the patches. 

(ii) Two-dimensional patchy habitat consisting of two patches 
with species having longitudinal dispersion in the second 
patch, while transverse dispersion in both the patches. 

(iii) Two-dimensional patchy habitat consisting of two patches 
with species having transverse dispersions different in 
the two patches. 
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6.3 EFFECTS OF MIGRATION IS' A OSS DIMEhSIQKAlj HABITAT : 

cto patches 

Consider the habitat 0 < x £ L to be divided into two 
linear patches by a line x = 1^ such that first patch is 
defined as 0 < x < and second patch as < x < Lg« 

It is assumed that the second patch is ecologically 
more favourable to both the species and hence the tendency of 
convective migration of the species in the first patch is 
more and dispersive migration may be negligible i.e. = 0, 
i = 1 , 2 . But in the second patch both convection and dispersive 
migration are dominant. It is also assumed that the convective 
velocities for the species are same in the respective patches 
i.e. ^ = Ug^ = U.., 3 = 1,2 while the dispersal coefficients 
are same (D^g = Egg = i n second patch. 

With these assumptions, the system ( 6 . 1 ) may be 
linearized about the equilibrium state by writing 

u i;j( x, t) = + v i3 (x,t), i, 3 = 1,2 


to get 


dv 

W 




av 


u. 

D 


ax 


ii = 


2 

I 

fc=l 




j = i, o < x < I 


(6.2-a) 
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+ v Hi _ | h „ 

dt Sx b + D :j = 2,Ij<s<i 


where 


( 6.2-b) 


af, 

b ik ~ ^au“^ 


k u i*“g 

lie system (6.2) may be associated with the following initial 
and boundary conditions 


y ij (x,0)=0; i,j = 1,2, x>0 

v 1;L (0,t) = H 1 > 0 

t > o 

v 21 (0,t) = P 1 > 0, 


v 12 ( Xft ) 0, 
v 22( x,t ^ °> 


t > C 


(6.3) 


(6.4) 


(6.5) 


as x - for semi-infinite habitat or as x - L, for finite 
habitat* 

The condition (6.4) implies that there exists a 
reservoir for both the species at x = 0. At the interface x = 1 
of the two interconnecting patches, the following matching 
conditions are prescribed : 
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Lt 

x — L- 

v l2 (x,t ) 

= Vu(Li,t) 


1 



(6.6) 

Lt 
x -* 

v 22 (x, t ) 

= ni CL i’ t) 



6 .3 .1 COI/rPESITIQi MODEL 


In the case of competition model, the interaction 
functions are given by, 

f lj = u lj (a l - °1 "zp 

(6.7) 

f 2 3 = U 2 3 ^ a 2 ” C 2 U 13^ 


where and u g ^; j = 1,2, denote the densities of two species 

in 3th patch. IDhe corresponding linearized system for the 
evolution of the two species are given by (6.2) with 


H* 

H* 

11 

O 

sm 

i = 1,2, 


12 - “ f “ 

- °lV°2 < ° 

(6.8) 


f,s > o 


'21 = = 

- C2 a l/ C l ^ ° 



(i) Semi- inf inite habitat 

Solving system (6.2-a) with (6.8) and using conditions 
(6*3) and (6.4), the distribution of the species in the first 
patch are obtained as follows i 

v u = Di MN H(t - V 


'21 


Y 21 


(S-) H(t - §-) 


1 


(6.9) 
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where 

= 1?^ cosh st - P 1 r sinh st 
V gl (t) = ccsh st - — sinh st 


( 6 . 10 ; 


s ^b 1 g.bgi , r - l/k-jj^Tbgg 


( 6 . 11 ) 


and H(t - ^-) is the unit step function. 

Similarly, solving the system (6.2-b) with conditions 
(6.3), (6.4), (6.6) corresponding to infinite habitat, the 
distributions of the species in the second patch by using the 
transformation similar to (*;*' : - , )» are obtained as 


v 1? = f <Mx-L 1 ,t-f)V 11 (t -1+ *1) H(: 


f?i) dl 
U 1 


v g2 = / <J> g (x-L 1 ,t-3!)V gl (t-f+ Tf i) H(2 _ dl 


( 6 . 12 ) 


wnere 


(x - ir-jty 


<f>.(x,t) = -jjg- . — ™ exp (- — — 
3 t^ Zinfr 


(6.15) 


and j=2. 


It may be noted from equation (6.9) and (6.12) that 

for TJ^ = 0, no species exists in any patch as expected from 

initial and boundary conditions. Prom (6,9) it can be seen th,t, 

X 1 

at any instant t = jj- , 0 the species in the region 

0 < x < x 1 have nonuniform spatial distributions while in the 
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region x > x^ the densities are zero. It may also he noted 
that the distribution of species in the first patch are unaffected 
by the presence of second patch. However, it can be seen from 
(6.12) that the populations in the second patch are very much 
dependent on the size of the first patch and the conditions 
prescribed at the boundary of the first patch. 


As t -* “ it can be seen from (6.9) and (6.10) that the 
limits of v 11 (x,t) and Vg 1 (x,t) approach to Y and 
V 2 i(x/Ui) respectively which are the steady state solutions of 
the system (6.2-a) showing that the equilibrium state is always 
asymptotically stable. 


From (6.12) it is noticed that as t - » the populations 


v 12 , v 22 would respectively tend to 


12 


(x) = V u (fli) exp ^^(x-l 1 ) - (g§ + s) 1/2 


22< x) = V 21^tjJ ) exp 2 § (x “ L 1 } ” ( 5§ + s)1/2 

provided 

S'! = P x r 


u; 


x-L. 


(6.14) 


v f 


(6.15) 


It can be verified that (6.14) is the steady state 
solution of the system (6.2-b) which would exist under (6.15) 
showing that the equilibrium state is asymptotically stable 
in the second patch also, but the convective and dispersive 
migration have no effect on stability. It is also noted here 
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that the stability condition (6.15) depends upon the boundary 
conditions prescribed at the first patch. 

Considering the particular case as 1^ -* 0 i.e. the 
single patch 0 < x < °° the equilibrium is asymptotically stable 
provided condition (6,15) is satisfied which is the same as 
obtained by Gopalsamy (1977). 

When, 1^ — °°, the habitat reduces to a single patch 
again and the equilibrium state is stable and the steady state 
distributions of the species are given by (6.9) showing non- 
uniform spatial patterns arising due to convective migration. 


(ii) Finite habitat 

In the case of finite habitat of length I, the solution 
in the first patch remains the same as obtained in (6.9) and 
the equilibrium state remains stable. However, in the second 
patch (1^ < x < l) the solution of the system (6.2-b), with 
initial condition (6.3) boundary condition (6.5) at x = 1 and 
the matching condition (6.6), by using the transformation ) 

is obtained as follows 


V lg (x,t) = 


T 22 (x >' t) = 


2D% 


°° pmit(x-L^) 

m — J 

00 l-mitte-L- )-i 

= J, m G 2i (I ‘i } sin L — J 


(i-i 1 ) 2 
2ln 


(1-1- ) 2 m=l 


(6 .16 ) 



I dtJ 
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where, 

Gudi) 


Xj Xj 

/ P ? (t-T, X -L 1 )Y ll (t-T+ £) HCT-^) df 


o 




Ti T) 

G 2 i(Li) = f I 2 (t-T,x-L 1 )T gl (t-I+ !T i; H(T - ^) df 


U„ 


^(^x) 


= exp 


Dm^-ji^t 




u? t tfj 

W ' UP + 


2D 




and j=2. 


( 6 . 1 ?) 


(6 . 18 ) 


It Is noticed from (6.16) and (6.18) that as t - 
Vl 2 (x,t), Y 22^ x,t ^ tend to 


2D7t 00 - mTc(x-L 1 ) 

Vf 2 (x) = — -* I m G 11 (L 1 ) sin ( — - _i-) 

1 (L-L-r m=l 1 1 L " L 1 


2Drc 


00 _ mn(x-L 1 ) 

v 22^ x) = 2 m G 21 (L 1 ) sin (-— -i- 

^ (I-L.) 2 m=l ^ 1 1 (L-L-l) 


(6.19) 


where 


Grn(ll) = exp 


i-U 


Li T. 2 2. u: 


jg(x- L 1 )].[jil((]i)(^ ^ + |f)-f V g 


( 6 . 20 ) 


21 


(1 1 ) = exp[ Joc-v] + 


3 't 1 - 1 ,-! 3 


2 2 ~ U? 

5 , - (T 2 -*^ + j$>‘ - - 


(6.21) 


_ 2 2 Uo 

+ 2 _ s > 0 




( 6 . 22 ) 


provided 
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It may be seen that these distributions are the steady 
state solution of the system (6.2-b) with conditions (6.5), 

(6.6) showing that the equilibrium state is asymptotically stable 
under the condition (6.22) and the stability increases as 
convective velocity or the length of first patch increases. It 
may also be seen from (6.22) that if convection is absent, the 
stability increases as dispersal coefficient increases. 


It may be noted again that even if the condition (6.22) 
is not satisfied, the limit of v^g(x,t) and v 22^ X} " fc ^ can ex i s ^ 
provided condition (6.15) is valid, giving 


v 12 (x) = 


2Dit 


U P 00 m L 1 i-m7i(x-I 1 )-r 

exp (gytx-ip)^ =- (_ 


v g2 (x) = 


H ) 1 


2Dtc 


U, 


21 


— T o ex P(of( x -Li)) 2 Z““ v 2i (fir)sin | — 

(1-I 1 ) 2 1 m=l J) 21 U 1 L ^1 


(6 .23 ) 


where 

2 2 tt 2 
_ m % D 

i 1 ~ o Tp + = 2 (6. 24 ) 

31 ( 1 - 1 , ± f 

which also satisfy the steady state system corresponding to 
(6.2-b) with boundary conditions (6.5) - (6.6) showing again the 
asymptotic stability of the equilibrium state under the 
condition (6.15), however, the stability does not depend upon 
convective and dispersive migration. 
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6.3.2 PR'BI PREDATOR MODEL 

Considering tne case of Yolterra type prey predator 
system the interaction functions may he written as 


f lj “ u lj ( a l - °1 “ap 


= ^ 0^2 + °2 U lp ! 3 = 1 ’ 2 


(6.25) 


where u^. and Ug.. denote the densities cf prey and predators 
in the jth patch. Then the equations governing the evolution 
are given by (6.2) with 


b ii 


12 


0 

-f 


f,g > 0 


(6.26) 


21 


(i) SEMI-INPIEITE HABITaT- 

The solution of the system (6.2) with (6.26) under 
conditions (6.3) - (6*6) in respective patches are given by 
(6.9) and (6.12) where ? 1:L ( t), Y 21 (t) are modified as follows i 

Y l:l (t)^ Y l2 (t) = N 1 cos (st) - P ± r sin (st) 

(6.27) 

N 1 

Y 21 (t) = Y 22^ = r sin + P 1 C0S 

In this case also as t -* ~ the population distributions 
in the first patch tend to nonuniform spatial steady state 



130 


distributions showing that the equilibrium state is asymptotically 
stable. 


In the second patch the limiting populations of two 
species as t - can be obtained from (6*12) and (6.27) as 
follows 


v lg (x) 


exp 


^-k 1 ( x -L 1 )J V 12 kg ( x— L ^ ) + 



(6.28) 


Vgg(x) = exp 




22 


lg(x-l 1 ) + 


s L. 

“U7 


a 


where 


kf = 


2 4 

u II 

2 ,1 / U 2 

— ~ + 2D ^ 

2T 16D< 


L-8D 2 


+ s 


2 s 1/2 "I 1 / 2 U 2 
J ~ otT 


2D 


(6.29) 


r u p i u ? 2 

2 + 2lF ^ p + s ) 

L 8D 2 ^ 16D 2 


2 s 1/2' 


V 


As the distributions (6*28) are also the steady state 
solution of the system (6.2-b) with (6.26), the equilibrium is 
asymptotically stable * 

Since (x-L^) decreases as 1^ increases therefore v^g, 
Vgg increases as 1 1 increases. Thus, it may be noted from 
equation (6.28) that larger the patch 1, higher the densities 
of preys and predators in the second patch. 

For the case = 0, U g = 0 the distributions (6.28) 
are the same as obtained by Hadelar (1974). 
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(li) FINITE HABITAT 


Consider the prey predator model gi^en by (6*26) in 
the case of finite habitat* The solutions In the first patch 
remains the same as (6.9) with (6.2?) while in the second patch 
the distribution of the species can be obtained as 


v i9 (x,t) - 2Dtc 


mn(x-l 1 ) 


L-et 


hi" 

1 - 


PTjir 00 pmis(x-L 1 )-r 


22 

where 


pm7i(x-Lj ) 
1 


t ii 1 L 1 

Gl 2 (ll) = / F 2 (t-l',x-L 1 ).V l2 (t-T+ |yi) u(T-^)df 


(6.30) 


(6.31) 


t 1- h 1 

Gg2^l^ = f ■^2^”"'^’ X "*‘^1^ *^" 22 ^"*^ u(T- u^*)dT 


and Fg(t,x) is the same as in (6=18). 

As discussed before , the equilibrium state is always 
asymptotically stable in the first patch. 

In the second patch the distributions in the limiting 
case as t -* «> can be obtained from (6.30) by noting from (6.18) 
that the coefficient of t is always negative, 
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2D% 


v (x) = ~ ±2 ~2 2 m G 12 (1 l^ sin ( "-T-T- -- 0 

(L-L. r m=l L<L 1 - Lj “ Jj - 


7 


(6 .32) 


v g8 ( X ) = 


2Dit 


00 _ m^x-l- ) 

<l*l’) 8 m=l “ G 2 2(LlJ S “ ( i- 1 ! J 


wnere 


G 12 (L 1 } ~ exp Ls§ (x - L l ) ] , L Yl2(t7 I ) * ( ^ir ) 


2_2 UN 1, 

& \ J? 7 T / JL 


2 + V 22^ ) J/ I 2» 

(6.33) 


G 22 (L 7 = eXp [s§ (x “ L l^ *[f Y 12(^ + ^~^y2 + S )Y 22^_^2 


which is the same as the steady state solution showing the 
asymptotic stability of the equilibrium state and the stability 
increases with convection. 


6.3.3 EFFECTS OF MIGRATION IN FINITE NUMBER OF PATCHES 

Let us consider that the habitat is divided into 
finitely many patches, say n. Then the linearized equations 
in this case are given by (6,2-a) where i = 1,2, j = l,2,...,n 
and jth patch corresponds to < x < L^ with 1 Q = 0. 

We discuss the following cases : 

(i) When only convective migration is present in all patches. 

(ii) When in the first (n-l) patches convective migration is 
present while in the last patch dispersive migration is 
also present. 
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In case (i) for competition model with interaction 
functions ( 6 - 8 ), the solution of equations ( 6 . 2 -a)with conditions 
(6.3), (6.5) and. the matching conditions at the interface x = L . - 

_L 

V i, = V i,^ L j-l^ 1 = lj2 ’ 3 = 2,...,n (o.34) 


can be obtained as 

Vq^Cxjt) = Y 1 ± (Ij(x)) H.. ( t,x) 

0 = 1 , 2 , .. • ,n ( 6 . 35 ) 

L . 1 < x < 1 . 

t (x,t) = y (T (x)) H (t,x) 3-1 3 


where 


x - L, . 

T ( X ) = ■ .. ,, ... , - 2 - 1 


for 3 = 1 


1 < D < n 


(6.36) 


3 u j k=l k 


and 


x - L . 1 

H.(t,x) = H(t y • ^~- ) for j = 1 


(6.37) 


c-1 . 


x- 1 . 


i L. L t 


= H(t $ 


J=i) S H(t.q^=i..£ 


D i=2 


3 k =2 “fc -1 


3 > 1. 


It can be seen from (6.35) that as t 
populations in the jth patch are 


v 13 (x) = V u (T 3 ( X )) 

= W T J< X)) 


°° the limiting 


(6.38) 
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which is the same as the steady state solutions of the 
corresponding system showing that the equilibrium state is 
asymptotically stable. 

Considering the case (ii) for the competition model 
again the solution in the jth patch for 1 < j < n-1 is the 
same as obtained in (£.35). In the nth patch of infinite extent 
the solution of (6.2-b) using Interaction coefficients (6.8) with 
conditions (6.3), (6.5) and the matching condition, 


v i,n-l ( L n.l } = T in 


(6 .39) 


is given by 


(6.40) 

t _ 

v 2n^ x ^ = ^ ^n^ x “’ dl n-l , ^”^^ , ^ r 21^”^ + ^n-l^n-l^^n-l^ ,d 'n-l^ dx 


where <l> n (x,t) is given by (6.13) when j = n and 
^n-l^n-l^ H(l,l n _i) are given in (6.35) - (6.37) respectively. 
As t - the population distribution in the nth patch can be 
obtained from (6.40) as follows : 


v 


r 


u. 


= exp | s#(x-I , )-(jS + ( 


In 


ll v ‘ n-l v n-1' 


L2D 


n-1' 


U 

_n 

•4D 


. 1/2 


x-L, 


n-1 


fu 


>] 


T 2n = T Sl( T n-l( I n-l )) ex P 



(6.41) 



provided (6.15) is satisfied. 
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It call be verified that the steady state solution of 
•the system (6.2-b) with conditions (6*5), (6.39) will exist 
provided (6.15) is satisfied and it is the same as (6.41) 
showing the asymptotic stability of equilibrium state which 
does not depend upon convective and dispersive migration. 

If the nth patch is finite then the distribution of 
species can be obtained from (6.2-b) and (6.8) as 


r .u-y 


v “ (x>t) = ^T17 ii" dncvh 3iD 


'n-lH 

IPJ 


(6.42) 


v p „(x,t) = 


“ r mit(x-L )-, 

^ m J 2n( L n-l } sil1 L 


' ~ , . T \2 


( L - L n-l) m: 


where 


O 

(6.43) 


l’ n (t,x) is given by (6.18) for j = n, 

and T lt (I!) and H 1 (1,L 1 ) are defined in (6.35) and (6.3?) 


respectively. 


As discussed before when t — °° the finite limit of 


(6.42) would exist provided 

~ 2 2 U? 

Dmn n 

ttz — 7S + SJ “ s > ° 


( L “ L n-l^ 


(6.44) 
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or condition (6.15) is satisfied giving' the steady state 
solutions in these cases respectively as follows 


V in (x) 


2£>7t 


C^n-P 


£ m J lP L n-P s “ 




1-1 


n-1 


(6.45) 


00 - mit(x-l 1 ) 

^P x ) = rff—v £ m ( -nrr- ) 


2Dn 


2n 

and 

7 iP x > = 




2m 


1 


ex P [sB^n-p] £ ll^n-l^n-P^ 

m 1 ^nl 


m 


c~VP 2 


pm*(:x-l ,}-i 

s “ L (t-t-pj 


V 2P X ) 


2Dlt 


t^-VP 2 


(6 .46 ) 

f— U Q -1 00 

exp [Jtx-vpj £ f-T 21 (Vi(VP) 


nl 

pnm(x-l 1 ) 
sin f ■ 


L H-P 


] 


where 


2 2 


U 2 

+ aS) 


? ln = ex » 55 [ I'-Vl 1 ] -fjlPVfVP^J* )2 ’ 45 

- f ^PVPVP^j^n 


^2p L n-P = 6X5 £ 55^ x-i n-1^3 * [^Sp^n-P^n-P^^ )2 + 35^ 

XI* - * l 

- S iPVPVPlK < 6 - 47) 
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D n and. D n ^ axe given by (6.21), (6.24) respectively wnen 1 = n 
which shows that the equilibrium state is asymptotically stable 
under the condition (6.44) or (6,15). In particular from 
(6.44) it is noted that stability increases as convective 
velocity or the length of first n-1 patches increases. Further 
in absence of convection , stability increases with dispersion. 

Finally, it can be remarked here that Case (i) is a 
one dimensional version of the problem studied in Chapter IV. 
However, Case (ii) may be regarded as a particular case of two 
subhabitats with first patch having variable convection 
approximating as (n-1) step functions while the second patch 
with constant convection and dispersion. 


6.4 EFFECTS OF MIGRATION IN TWO DIMENSIONAL PATCH! HABITAT : 
TWO PATCHES 


Consider the following particular cas® of two competing 
species in a two dimensional heterogeneous habitat consisting of 
only two patches when 

1. Convective velocity of both the species in the two patches 
are equal, i.e. 


U. . = U. 
ID 3 


V. . = V . 
3-D D 


i = 1,2? j = 1,2 


longitudinal dispersion in first patch is zero for both 
the species i.e. 


D 


ilx 


0, i = 1,2 . 


2 
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3. Dispersion for both the species are equal in two patches 

i.e. 


D i2x - D 1 ’ I) i2y " D 2 


D ily " 33 

Under these assumptions the governing linearizeu system 
(6.1) for competing species can be written as 


av. 1 av.. av 2 

FT 1 + u l + v l ^ b ik + D “2 


-2 

a v 


il 


(6.48-a) 


3 v i2 

FF“" 


av. P av. p 2 a 2 v, p a 2 v ? 

+ TJn + V c = 2 b v 2 + D — ^ + D — g*- 

3e=l 115: ^ - 1 dx 2 ^ ay 


2 ax T '2 9y 

(6.48-b) 

with interaction functions (6.7). 

The linearized system (6.48) is solved under following 
conditions 


(x,y ,0) = 0 ; i, 3 = 1,2 

(6.49) 

'lj(x,0,t) = IL 5 v g; .(x,0,t) = P^ 3 = 1,2 

(6.50) 

i; j(x,y,t) - 0 as y - B ; i, j = 1,2 

(6.51) 

i i(0,y,t ) = 0 ; i = 1,2? 3 = 1 

X J 

(6.52) 

1i (x,y,t) -* 0 as x -* 1 , i = 1,2 ; 3 = 2 

X J 


It v i2 (x,y,t) = v il (l 1 ,y,t) ; i = 1,2 

(6.53) 


x -* L 
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Since in this case the situation in the first patch 
is same as discussed in Chapter V, the solution of (6.48-a) 
with conditions (6.49) - (6.52) are same as (5.18) - (5.22) 
and the condition for asymptotic stability in the first patch 
is given as 




s > 0 


or 


N i = p i r 


Again the solution of the system (6.48-b), (6.53) in 
the second patch can be obtained as, 


v i2 (x,y,t) = J i2 (x,y,t) + J i2 (x,y,t) (6.56) 

where 

4TtDp 00 . 

J.„ = - — - I 2 n . G. (t,m,n) e(m) 

1 LB 2 m=l n=l 


m7i;(x-L 1 ) 


QJW 

B 


(6.57) 


J i2 


4nD 1 °° °° 

2 2m G?(t,m,n) 
I 2 B m=l n=l 


sin 


m7t(x-l 1 ) 


G*(t ,m,n) 


1 

/ 

-U, 


si n£K 


(6.58) 


/ V 12 (I) e X p[ 5 |-( X -L 1 ) + 3 §- y -F 2 (m,n)l_ 


dS 


(6.59) 
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G?(t,m,n) 


2tiD 

B 2 


n 

2 i: I(lc,n) . 
k=l 


t t-I r U 9 

/ / ( X -L 


} + !g_ . _» 
l} *=2 


t*““ ij-j •“ 


P 1 (]£)T 1 -B 2 (m,n)T^j 

(6.60) 


dl^dl 


M*) 


l Y 1 


B* 


535 


P 2 (m,n) = — sg-i + 


2 2 _ 2 2 -^ tt 2 —2 
m it n tt D 2 U g V g 


2 + 4ET + 4^7 


(6.61-a) 
(6. 61- Id ) 


0(m) = § 

L 


mn 


«p (- sfifj L^t- + 4 ^J 


, Vi-irA* 4 d- 1 

(“ 9PT“ — 5“ + — 5T 


I(k,n) = J' exp 
o 


( 


Zi Y s_ 

2D " 2D, 


*] 


sin Sffi. ain SS 


B 


(6.62 ) 


sin ^ dy (6.63) 


1 = 1 - 1 1 

and H(t) is the unit step function. 

In the second patch, it can he noted from (6.56) - (6.60) 
that the distribution of species is very much dependent upon 
the distribution of species in the first patch. As, t - ■*>, 
finite limits of v^ 2 , i = 1,2 would exist provided J^ 2 and J^ 2 
are finite as t - °°. Let us explore the possibilities under 
which this may happen. 


It can be seen from (6.57) that, as t - the limits of 

J i2» 


i = 1,2 exist, provided 
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L 2 B 2 


U 
£ 

+ TrT" + 


T 2 

J 2 


w: + w: - s > 0 


(6.64) 


or 




> y» 
2 x 


(6.65) 


and would be given by 

4-n3 0 


Lt J 

“j) — > oo 


IP = =-p~ exp 
1C L3^ 


u 0 v 2 ■ 

1 ^ x “ L i^ + 2^y 


_2 

23 


oo oo 


2 2 


ne (m ) 


m=l n=l Fg(m,n) 


" fP 2j * 


sm 


■ nm(x-L^) 


. sin SgSL 


( 6.66— a) 


1 t Jq Q — 
t - « ^ LB 


4-nDp r U ? V„ - T 

* exp ( x ~ L i) + 23^ y J * 


oo oo 


2 2 r 


ne (m) 


m=l n=l Bg(m,n) 

irm^x-L^) 


[p 2 E 2 (m,n) - gHg] . 


sin 


sin 


(6.66-b) 


where 


Bp(m,n) = p?(m,n) - s 2 


(6.67) 


or 


Lt J 


12 


47uD f 


r ir P v 2 *1 

• e ® L337 (*- L P + sr y_ 


ne(m)Up 

2 — 


oo oo 
2 

m=l n=l B 22 (m,n) 


dot(x-L 1 ) 

sin sin 


(6 .68— a) 


nity 

B 
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Lt J 


22 


4tlD, 


IB* 


exp 


u2«l 


(::-L 1 ) + 


V 2 ' 


00 00 n0 (m) P 

2 2 


m=l n=l F 2g (m,n) 


init(x— L 1 ) ___ 

-sin — — i- sin 

T 


(6.68-b) 


where , 


P 22 (m,n) = F g (m,n) + s 


(6.69) 


under conditions (6.64) or (6.65) respectively. 

However, if longitudinal convective velocity is zero in 
the first patch then the finite limits of J ig as t - ® would 
exist provided 

(i) (6.54) and (6. 64) are satisfied, or 

(ii) (6.55) is satisfied 


and the corresponding limits are given by 

8 AlL p TJ 0 V 2 1 

1 Ji2 = ' exp (x_1 i ) + ^ 7 j • 

(6.70-a) 

0 ° co ml. (m,n) m7i(x-l 1 ) 

2 2 — i sin ™ — — sin ^ 

m=l n=l P 9 (m,n) L 

where 

1.(11,!!)= s h^he2r F (jc) (s 1 r„(m, 0 ) - f Pji 

± te=l I' 1 (k) >- 

- f {P 1 l' 2 (m,n) - g Nf } J 


(6.70-b) 
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I 2 (m,n) = I 


k I(k,n) 


k=l (k) 


k 1 (k) {P 1 P 2 (m,n) - g N 1 } 


g { U 1 k 2 (m,n) - f P 1 > 


(6.70-c) 


i,(k) = pf(k) - s" 


(6,70-d) 


or 




87t 2 DD 1 | 

r Up(x~L 1 ) 

Y 2 1 

Lt 

— *■ oo 

J 12 " 

exp 

ii X) 1 

L-VB7 

+ 



00 00 m ^(n)!-. 

ma(x-L 1 ) 



2 r 

— — - Sin 

= — — sin 



nt=l rt=l PggCm^n) 

L 


_ 

p. 



Lt 

J 22 = 

IT ... Lt J 12 


-*■ oo 


X T> -> oo 




mm 

00 k I(k,n) 



I (n) 

= 2 — 





k=l f' 1 (k) 

+ s 



(6.71-a) 


(6.71-b) 


(6.71-c) 


When longitudinal convective velocity is nonzero in the 
first patch then the finite limit of Jig as t -* » would exist 
again provided (6.64) or (6.55) is satisfied and the limits 
can be obtained as follows : 


It J 
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8h 2 DD 1 

=TT- exp 
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r u 2 , T w v i v + y 2 v -i 
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(6.72-a) 


^ . exp (6<re _ b) 
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nt=l n=l Pg(.m, 


* sin 


m7r(x-l 1 ) 


sin 


nrc.v 
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I 3 (m,n) = 2 k I(k,n ) £F g (m,n) & 1 (L 1 ,k) - f G g (L 1 ,k)]] 


( 6 . 72 - 0 ) 

oo 

I 4 (m,n) = ^2 k I(k,n) £F g (m,n) ^(L^n) - g G 1 (L 1 ,k)^| 

(6.72-d) 

mm 

G i (L 1 ,k) = f V il (2/U 1 )exp^ y - P ± (k) Z/U^JdZ (6.72-e) 


1 c 


under (6.64) and 

87t*TUL 


t J 12 = E 5 ? 


[>- 


i 


Ur 


e xp (- ^ J .exp [_^(x-L 1 ) + ^ y 

(6.73-a) 


OO oo 

2 2 


m=l n=l I , 22 (m,n) 

P, 


eJJjlL Si sin 4 sin «i 


Lt J 22 = F7 . Lt oo J 12 

tz CO X) —► CO 


(6.73-b ) 


when (6.55) is satisfied. 

Keeping in view of finite limits of J^ 2 and <K g under 
different conditions as pointed out above, it can be noted 
from (6.56) that finite limits of v 4g (x,y,t) i = 1,2, as t — =° 
would exist under one of the following alternatives, for = 0, 


Al i (6.54) and (6.64) are satisfied 

A2 : (6.55) and (6.64) are satisfied 

A3 : (6.55) and (6.65) are satisfied. 


Similarly the finite limits of v ig (x,y,t) as t - °° 
under one of the following alternatives, for ^ 0 


would exist 

t 
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A3 i (6.55) and (6.65) are satisfied 
A4 : (6.64) is satisfied. 


Wow, it can be verified that under above alternatives Al - A3 
and A3, A4 the limiting population v^ 2 (x,y,t) as t - °° are the 
same as the steady state distributions v? 2 ( x >y) °^. corrss “ 
ponding system given by 


v?g(x,y) = 


for U 1 = 0 and f 0 respectively, showing that the equilibrium 
state is stable in the second patch. 

However, if none of these alternatives is satisfied, the 
equilibrium state in the second patch is unstable. 

When longitudinal convective velocity is zero in the 
first patch and (6.54), (6.55) are not satisfied, it can be 
noted from Al - A3 that the equilibrium state which is unstable 
in the first patch remains so in the second patch also. It is 
also seen that even when the equilibrium state is stable in 
the first patch, it may not be so in the second patch. 

When 0, the situation is fairly conplex and whether 

the stability of the equilibrium state in the second patch 
depends on the first patch or not would be decided by the 
alternatives (A3 or A4) satisfied by the system. In particular, 
for U 1 0 the equilibrium state which is always stable in the 
first patch may remain so in the second patch also (see A4). 
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It is also noted that for sufficiently small and nonzero 
■values of dispersion coefficients and nonzero convective 
migration the condition (6,54), (6,55) can always be made to 
satisfy and hence the equilibrium state is always stable In 
both the patches. 

From the above discussion, it may be concluded that for 
competition model, the patchiness of the habitat may cause 
instability of the equilibrium state wnicn may be stable 
otherwise. 

6.5 EFFECTS OF PATCHINESS DUE TO TRANSVERSE DISPERSION IN A 
TWO DIMENSIONAL HABITAT j TWO PATCHES 


In the following , another particular case of competing 
species in a two dimensional finite habitat is considered when 

1. Transverse convection for the two species in both the 

patches is zero V . . =0, i,j = 1,2. 

d 

2, Longitudinal convection and dispersion for the two species 
are equal in the two patches i.e. 

= u * na = 1 > 2 

; i,a = 1.2 

Under these assumptions >the governing linearized system 
for competing species can be written as 
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with b^. defined in (6.8). She linearized system is associated 
with zero boundary conditions and nonzero initial conditions 
and matching conditions (4.4,4.5-a)-(4*7)as discussed in 
Chapter IV. 

Using the following transformation 


V lD 





r cosh st 
-sinh st 


sinh st 
-1/r cosh si 



“7 “i 


1 Z 9 . 

_ 2 3J 


,U IT , v 

exp (S3? x - IE ^ 


(6.75) 


the system (6.74) transforms to (4,11) with zero boundary 
conditions and matching conditions (4. 13,4. 14)Hence the final 
solution obtained in this case will be given by (4.23) multiplied 


by exp ( 


U 

2D 


x - 


i r 

4D 


t ). 


Thus the condition for stability for the 


equilibrium state is 



.T * 2 T 
+ 1 


2l 

2 


- s > 0 


(6.76) 


where X^, Xg are given by equations (4.26-a) and (4.26-b)showing 
that longitudinal convection further stabilizes the equilibrium 
st ate . 


In the case of prey predator model also the effects 
of convective and dispersive migrations due to finite number 
of patches in the habitat can be investigated and similar 
results as pointed out in the case of two patches may be found. 
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6. ,6 CONCLUSIONS 

In this chapter the linear stability of migrating speci 
in a patchy habitat has been discussed by considering various 
combinations of convection and dispersion in different patches. 

In the case of one dimensional habitat with dispersion 
only in the last patch, it has been shown for competition model 
that an inherently unstable equilibrium state becomes stable in 
all patches in absence of dispersive migration. However, in the 
case of semi-infinite habitat the migration has no effect on the 
stability of the equilibrium state even in the presence of 
dispersion in the last patch. But when this patch is finite, 
a general condition for stability involving convective velocity 
and dispersion coefficient has been obtained. 

In the case of two dimensional habitat consisting of 
only two patches, the linear stability of interacting system has 
been discussed with dispersion only in the last patch while 
transverse dispersion in both the patches. In the case of 
competition model, it has been shown that the stability of 
equilibrium state in the second patch is dependent on the 
conditions of its stability in the first patch and even if the 
equilibrium state is stable in the first patch, it may not be 
so in the second patch. But if the equilibrium state is 
unstable in the first patch, it would always be unstable in 
the second patch. | 
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However, in the presence of equal longitudinal 
convection and dispersion of the species in all the patches, it 
has been shown that the effect of longitudinal convection is 
stabilizing in absence of transverse convection. 

The case of prey predator model, in all these cases 
has also been discussed and it has been shown that the equilibrium 
state is stable with convection and dispersion in patchy habitats. 
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CHAPTER VII 


EFFECTS OF SELF AID CROSS DISPERSION OH LINEAR 
AND NON LINEAR STABILITY 


7.1 IN TRODUC II ON 


In previous chapters, we have studied the linear stability 
of equilibrium state of two interacting and migrating species 
systems in finite and infinite habitats by considering variable 
dispersion coefficient. However, it may be noted that, 
biologically, it is more important to study the stability of 
nonlinear systems even if it cannot be solved by analytical 
methods. Such nonlinear stability analysis with and without 
dispersion has been investigated by using Liapunov second 
method in few cases (Gatto et , al., 1977, Jacob Jome and 
Shlomo Carmi, 197 7; Gopalsamy et. al., 1980). 

Keeping the above in view, in this chapter, we study the 
linear and nonlinear stability of two interacting and migrating 
species system by considering density dependent self and cross 
dispersion in a finite two-dimensional habitat. 

7.2 BASIC EQUATIONS 

A general model for the evolution of two interacting and 
migrating species in a habitat, 0 < x <. L, 0 <_ y < B, in absence 
of convection, can be written as [see Chapter II[]» 
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3u 


eu. 


vr = f f u ’ u 2 ) + h D^oVsr + dlx s 


Siu dUp_ 

2+ |y2 ) iy^ a ^^7‘ +d iy~8y-* 

(7.1) 

su g au.^ 

au 3u 2 a yE > 2y (u 2 ) 8y _+a 2y aF- 

at 2 = f g( u i> u g) + & 2 i C u 2 ) ax 2x 


U V ^ U.A — * Cn . 

) the self dispersion 

where u ± , i=l,2, the density; ( B ix’ JJ coefficient and 

■\ fVip cross disp^ r 

coefficient and (d ,d J th SDeC i e s. It 

i - function of the 1 sp 

f . (u-,u ) is tile interaction native functions 

12 4-1.2 are positive 

T > » 1 L m As 


fuTiCbXO^ UJ ~ 

f . (u-,u ) is tile interaction native functions 

12 . i= l,2 are P 0S1WVB 

has been assumed that D ix > D iy’ . which may be 

• *1 £ are cons 

of u ,u while d^, d^; 1 ’ oI interaction 

upon the iyp o 

positive or negative depending v 
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sen the species. . . r \ 

rontritial e 4 uiliW iu “ <*»* 

The positive uniform nont 

v+ , ned by solt^S 

for the system (7.1) can be obtaxne 

>1 o 


mw 01 1=1,2 

using the transformation, 

1 = 1,2 

u. = u. + v 4 * 
i l i 

the system (7.1) can be written as 
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0v, , 

8v i /- - x 3 fi> (u-d+v ax" + Sx 
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Linearizing (7.2) we get 


2 a 2 v,. _ aV, a 2 v a 2 v, 

^ + IT" + TT + ^lx — 2~ + ^°v — 9* 

c=i x x 0x^ ~ y 1 sy* lx ax 2 ay 2 

(7.3) 


1 ^ 1 b 2k V L + -W V d p 


a 2 a 2 a 2 

3 v ? a v- a v 1 

+ W ^S" + d 2x ^5- + d 2y rr 


where, 

af, 

hk = <au^ 3 i’ u z > i>k = 1 > 2 

It can be noted here that in the linearized system (7.3), 
the self dispersion coefficients are now functions of equilibrium 
state and' would be constants for fixed (u^,Ug). 

The following initial boundary conditions are associated 
with the system (7.2) or (7.3) 


v j_(x,y ,0) = F i (x,y) 

\(0,y,t) = 0 

v j_(i »y , t) = o 

v j_(x,0,t) = o 

v i (x,B,t) = 0 


(7.4) 


(7.5) 


In the following we discuss the linear and nonlinear 
stability of equilibrium state in cases of competition and 
prey ' predator models. 
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7.3 COMPETITION MODEL 


In case of competing species the interaction functions 
are written as 


f l u l^ a l“ b i u l" c i a 2^ 
fg = Ug(ag- b g u g - c 2 U 1^ 


(7.6) 


The positive equilibrium state of the system (7.2), (7.3) 
with (7.6) can be obtained as 

- a l b 2 ~~ C l a 2 . - a 2 b l~ C 2 a l 

U 1 ~l b "2" C i C 9 ’ b 1^2“ C l C 2 


provided, 

b l v a l v C 1 
c 2 a 2 ^2 

or, 

°2 a 2 ^2 

The corresponding linear system is given by (7.3) with 
b. 

b. 


(7.7-a) 


(7.7-b ) 


ii = -Vi < 

0 

. h 

J 12 

21 = "°2 U 2 < 

0 

. h 

* 22 


■ C 1 U 1 < 0 


(7.8) 


-b 2 u g < 0 


7.3.1 DISCUSSION OE LINEAR STABILITY 


To discuss the linear stability of the system (7,3) about 
the equilibrium state (uL,u.g), let us choose a positive definite 
function, E, in the first quadrant, as 
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„ 1 B r h , 2 2 N , , 

B - 9 J J (v. + cv ) dx dy 
^ o o - L * 


(?. 9 ) 


whore c is a positive constant tc be determined. 
From (7.3) and (7.9), aE/dt is obtained as 


ff = J J L b ll V l + cb 22 T 2 + ( b 12 +cb 21 } v i v 2 ^ dxdy 


B L 


o o 


B L 


*2 
3 v 


a 2 
3 v. 


.2 

3 Y r 


a 2 
3 Yr 


1 . „ u V 1 , ^ _ v a 2 . _ „ u 2 


+ { i + h/i -r + ^ ^ + o3 V* ^ 

3 2 y ? 9 2 y 2 3 2 v 

4- d. ^ + d v. — + cd 0 v — rr^ (7*10) 

lx 1 ax 8 ly 1 ay 8 2x 2 ax 




+ cd 0 y, 


2y 2 ay 2 


d dz '-iy 


Using zero reservoir boundary conditions (7.5) around the 
habitat the second term on the R.H.S. of (7.10) can be simplified 
giving 

B L 

II = / / D>u t ? + ob 2 2 v I + ( b 12 +cb 2 f V Z 2^ 


B L 3v. 2 av. 2 , 3v p 2 3v p 2 

i / LP^aF") + D ly^3y“^ + 0 {j} 2x^3x"^ + D 2y^3y“^ 


(7.11) 


3 v 1 3v Q 3v 0 3v 1 

+ ( d i x +cd 2 X )(ax" ^6x“ ^ + ^ d lv +cd 2v ^3v”^3y~’^ dxdy 


Denoting 


P = b^v* + cb 22 v 2 + (b 12 +cb 2l ) rf z 


(7.12) 
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av 1 2 9v, 2 __ 3v 9 2 8Vo 2 

Q “ ^lx^Fx”'* + ^ly^dy”) + C ^2x^8x“^ + 3 2y 


3v 3v 6v Qv^ 

+ (d 1 + cd )(_— 1)(— -S.) + (d + cd,. Xtt'—Kt — 

v lx 2x 3x v 3x ly 2y 3y 6y 

we get from (7.11), 


(7.13) 


j-p 3 L B L 

ff = / / Pdxdy - / / Qdxdy (7.14) 

OO 0 0 


In absence of dispersion, it can be noted from (7.12) 
that E is a Liapunov function provided the quadratic P 
is negative definite, the condition for which is obtained as 

b ll b 22" b 12 b 2l > ^ b 12“ Cb 21^ 2 - 0 (7.15) 

choosing c = c^/cg we get the necessary and sufficient 
condition for stability of the equilibrium state without 
dispersion, after using (7.8) as 

b l b 2 ” c l c 2 ^ 0 (7.16) 

In the case of dispersion, when the condition (7.16) is 
satisfied, it can be seen from (7.14) that for E to be a 
Liapunov function, the other quadratic Q should also be positive 
definite, the conditions for which are obtained as 

4c DlA* > + 


4c D ly D 2y > ^ d ly + cd 2y^ 


and 


(7.17) 



156 


which can be rewritten as 

^lAx-'hx'W > (d lx" cd 2x )2 2 0 

and (7.18) 

4c(Di y D 2y ~ d i y d 2y ) > (d^-cdgy) 10 > 0 

Since D. , D. ; i=l,2 are positive .these conditions wo^id" 'hyty* 
xx iy 

be satisfied provided one of the following holds good. 


i ♦ 




d lx d 2x> " l "lx "2x1 


d 


xi, 


■iy 

d 2yl 

< 

l D ly V 


d 

< 

o 

‘lx 

2x 




d 

< 

0 

'iy 

2y 



in 

this 

case the cross 


(7.19-a) 


(7.19-b) 


of both the species are positive, the condition for linear 
stability of an otherwise stable equilibrium state is given 
by (7,19-a). However, in the absence of cross dispersion 
the otherwise stable equilibrium state remains so with seif 
dispersion. 


In the case when the quadratic P is not negative definite 
i.e. the condition (7.15) is not satisfied then with positive 

definite quadratic Q (i.e. the condition (7,19) is satisfied), 

% 

there exist a possibility for which dE/dt may become negative 
definite. In the following we find the conditions for which 
E becomes a Liapunov function in such a case. 
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As pointed out earlier the cross dispersion coefficients 
for competing species are positive, the expression for dE/dt 
can then be further, simplified to give 

B L 


dE __ 
dt 


/ J jP-CL 2 dxdy 


(7.20) 


o o 


where 


GL 9v 1 2 
x \ / r 


G 9v 0 2 G . 
Xu 2\ , X r a 


Q 1 " (I) lx" 2 ^ax ^ + ^ D 2 x C “ 2 ^8x ^ + 2 f 3x ( v l +v 2^ } 

G.. avu 2 G 3 Vq 2 G ^ 2 

+ ( D ly- + ( D 2y 0 - ■S £) (— > + 2* <ay (v l +y 2 )} 

(7.21) 


(7.22) 


G = 

X 

d lx + 

Cd 2x 

and 

G y = 

d iy + 

cd 2y 

considering the 

case 


B i x > 

G x ; 

cD 2x 

D iy > 

a 

y ? 

cD 2y 


(7.23) 

Gy 

for which the conditions (7.18) and so (7.19) also are consistent 
and using the Poincare’s inequality (Penn, 1975 ), (7.21) is 
simplified as follows. 

ft i { { C ( rt 22-°2 ) v 2 + (b 12 +0b 2l- S) 7 l 7 8^ 4xay 

'(7.24) 

where 

G = p 1 + cpg (7 *25) 
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ix ^2 


+ d. 


71 


xy -p.2 


a . 
l 




2 

D. ii_ 

xy b 2 


(7.26) 


It can be noted from (7.24) that the function E is a 
Liapunov function provided 

b ll b 22" b 12 b 2l'' < ‘ b ll CT 2 +h 22 CI l' 1 + ( b 12 p 2 +b 2l P l^ 

f 2 

+ (V2 " 4c. ) > 0 (7.27) 

It can be seen from (7.25) that 

G 2 > 4cP 1 Pg (7.28) 

and thus the condition for stability (7.27) takes the 
following form: 

b ll b 22“ b 12 b 2l“^ b li a 2 +b 22°l' ) + ^ b 12 P 2 +b 21 P l^ 

+ ( a 1 a 2”P 1 P 2 ) > 0 (7.29) 

under (7 .23) 

It can be remarked here that the condition (7.29) is 
necessary and sufficient for the stability of the equilibrium 
state under (7.23). Thus the equilibrium state which is stable 
or unstable with or without self dispersion would become stable 
with cross dispersion under (7.23) provided the condition (7.29) 
is satisfied. 
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7.3.2 DISCUSSION OF NONLINEAR STABILITY 


To discuss the nonlinear stability of the equilibrium 
state of the system (7.2), consider the following positive 
definite function in the domain u^ = u^+ v^ > 0, Ug= Ug+Vg > 0 


V = / / [v-j-i^ In (1+ -~) + c (Vg-Ug In (1+ -~)} Jdxdy 

oo u x 5g 

(7.30) 

where c is a positive constant defined as before. 

In the following we find the condition for the function V 
to be a Liapunov function providing the sufficient conditions 
for the nonlinear stability of the equilibrium state for the 
system (7.2) with (7.8). 


Let us first prescribe the function for the self dispersion 
coefficients (D- ,D. ) i=l,2 as 


D ix 


D . = c . . u. 

lx n i 


D . = D . + © . - u . 

iy ly ii i 


; i=l,2 


(7.31) 


where c..,e.. are positive constants, 
n’ ii v 


The time derivative of, V, can be obtained as 


dV 


B L 


= / 
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/ G — zr 

o v 1 +u 1 


8v 

Ft 


i + 


cv r 


9v r 


V u 2 


w -i dxay 


(7.32 ) 


Substitution of (7.2) in (7.32) would give 



160 


dV 

dt 


= - Sf (P + Q) dxdy 


(7.33) 


where 
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(7.34) 
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Sv 2 2 
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‘2“2x , / uv lw uv 2x . r Tl^l.y . w “2~2y 
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sv. 3v Q 

(gy-)(ojr‘) 


(7.35) 


In the absence of dispersion, it can be noted from (7.34) 
that V is a Liapunov function provided the condition (7.16) 
is satisfied giving the condition for nonlinear stability of 
the equilibrium state, which is the same as obtained for the 
linear stability , 


In the case with dispersion, it can be noted that if the 
quadrative Q is also positive definite then V is a Liapunov 
function (with self and cross dispersion) the conditions for 
which are obtained as 

j- 5 HX 5 2x~ d lx d 2x + °11 C 22 + °11 5 2 x 

1 2 (v 1 +u 1 ) 2 (v 2 +u 2 ) 2 (v+u 1 )(v 2 +u 2 ) (u^Kvg+Ug) 2 

+ — > Q. 

(Vf+Ug ) (^2+^2^ 


(7.36-a) 
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and _ 

u u 13 D iy p 2y~ d i.y d 2.y + e ii e gg + _ e n I) 2y 

1 2 (v 1 +U 1 ) 2 (v 2 +U 2 ) 2 (v^iL^Vg+Ug) (v.j+u^)( Vg+Ug ) 2 


+ 3 > 0 (7*36-a) 

(vi+Ui)^(v 2 +a 2 ) 

which are satisfied provided, 


^ly^2y“ d ly d 2y > 0 


(7.36-b) 


As these are also the conditions for linear stability of the 
equilibrium state, it is concluded that the sufficient condition 
for nonlinear stability of equilibrium state with dispersion is 
that it is linear stable with dispersion. 

If the condition (7.5) is not satisfied, then it can be 
noted from (7*33) that under the condition (7.36), there exists 
a possibility for which the function Y becomes a Liapunov 
function in a subregion of the positive quadrant. For this 
let us choose positive constants U o ,Y o such that 

^o > U 1 5 ? o > u 2 

In the region A = { (u^Ug): u x < 1°; Ug < Y q y . The expression 
for dV/dt (7.33) simplifies to 
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where 
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As discussed before (d^jd^y) 5 i=l,2 are both positive 
in this case becomes 
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where 
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(7.40) 


considering the case 
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(7.41) 
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for which the conditions (7.36) are consistent, we get after 
using Poincare's inequality 


av 

at 


B 1 uL u 

-f f c ( t) i + ” <n, + 


o o 
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qr x * 

o o 

Thus, the condition for stability of the equilibrium state can 
be obtained as follows 

- W 1 - a p W 9 

b l b 2~ C l C 2 +b 2 U l^ + TT^ + u 2 b l ^ + H 

0 o 

_ 0-i w. 0p Wp CpU-i C 1 U-p . 

+ U 1 S 2 ( ~2 + TT^ + vr } " ( "^“ P 1 + “^T P 2 } 

u o 0 v o 0 v o v o 


U 1 U 2 P 1 P 2 
Ij2 V ^ 

o o 


> 0 


(7.44) 


Thus, it can be concluded that (7.44) is the sufficient 
condition for nonlinear stability of the equilibrium state with 
dispersion in a region A under (7.41) even if (7.5) is not 
satisfied. As (7.44) determines the region A, it may be pointed 
out that the region of stability increases with increasing 
self dispersion (constant) and decreasing cross dispersion. 
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7 . 4 PREY-PREDATOR MODEL 


In the case of prey predator model, the Vol terra type of 
interaction functions are given by 


f-^(u^,Ug) — ( a.^— b .^u^— c .^Ug ) u 1 

f 2 ( VV = ( - a S- b 2V°2 U l )a 2 


(7.45) 


which gives the positive equilibrium state as follows; 


a 2°l + a 
~z + b 



cl-| 

C^^b^bJ 


provided a ± c 2~ a 2 1 °l > ® * The linearized system can be 

given by (7.3) with 


'll 

= -b^u^ 

< 

0 

b 12 “ 

_C i a l < 

0 

21 

= + °2 5 2 

< 

0 

b 22 = 

“ b 2 U 2 > 

0 


(7.46) 


To see the linear stability of the equilibrium state, the 
same positive definite function E is considered from which, 
proceeding as before, it can be noted that -E given by (7*9) 
is a Liapunov function in the absence of dispersion, as the 
condition (7.16) is always satisfied. Thus, in the case of 
prey predator model the equilibrium state is nonlinearly stable 
without dispersion, further, in presence of self and cross 
dispersion the quadratic Q remains the same as before and the 
condition (7.19) for positive definiteness of Q is always 
satisfied as the cross dispersion coefficients for the prey 



165 


predator species are of opposite signs (Kerner, 1959 ), Thus 
it is concluded that prey predator model is linearly stable 
with self and cross dispersion. 

fo investigate the nonlinear stability of the equilibrium 
state, choosing the same Liapunov function as (7.30) and 
proceeding, as in the case of competition model, it can be 
concluded that prey predator model is nonlinearly stable with 
or without dispersion. 

7.5 CONCLUSIONS 

In this chapter, the linear and nonlinear stability of 
uniform positive equilibrium state for both competition and 
prey predator models have been investigated by the Liapunov 
second method. 

The following results have been obtained in the case of 
competition model. 

1. The linearly stable equilibrium state remains stable with 
self dispersion. The necessary and sufficient conditions for 
the stable equilibrium state to remain stable with cross 
dispersion also has been obtained. 

2. The equilibrium state which is otherwise linearly unstable 
may become stable with dispersion (self and cross) under 
condition (7.29) for the case (7,23). 

3. In absence of dispersion the linearly stable equilibrium 
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state remains nonlinearly stable. 

4. In presence of dispersion it has been shown that the 
equilibrium state is nonlinearly stable in positive quadrant 
( u i>Ug) provided it is linearly stable with and without 
dispersion* 

5« However, if the equilibrium state is linearly unstable 
without dispersion which becomes linearly stable with 
dispersion (self and cross) in positive quadrant, then it 
is possible to determine a subregion where the equilibrium 
state is nonlinearly stable and the region of stability 
increases with increasing self dispersion and decreasing 
cross dispersion* 
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CHIP TER VIII 


EFFECTS OF PATCHINESS OB LIB EAR ABB IOHLIHEAR STABILITY 


8.1 IBTRODUCTIOB 

In the previous chapter, the linear and nonlinear stability 
of the equilibrium state has Deen studied with density dependent 
dispersal using Liapunov second method. As pointed out in 
Chapter IV, the dispersal coefficient may also vary with space 
and such effect on stability can be studied by considering the 
patchiness of the habitat. 

Keeping this in view, in this chapter, we study the linear 
and nonlinear stability of the equilibrium state of competition 
and prey predator models in a two dimensional finite patchy 
habitat. 


8.2 BASIC EQUATIONS 


Assuming the habitat to be consisting of finitely many 

th 

patches (say B ) such that the 3 patch is defined as 
L. < x < L.; 3=1,2, ... ,N with L o =0 and L^= L, the system 
governing the evolution of two interacting, dispersing species 
in 3 th patch can be written as 



dZu ii 8 %.i 

f - J + Us* + D i jy -CT 1 


ID 


ax 


ay 


(8.1) 


i = l,2j 3 = 1,2 ,... , 


M. 
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where iu.., (Ih^, D^^), f^.. denote the density, dispersion 
coefficient and the interaction function of i"^ species in 
patch. As before the interaction coefficients of the species 
have been assumed to be unaffected by the patchiness of the 
habitat and then the uniform positive equilibrium state 
(Ui»5 ? ) is obtained by solving 

f ij < 5 1>V = 0 

Writing 

u. . = u.+ v. . 

ID 1 13 

in (8,1 ) we get the non-linear form as 


av. . 


a 2 

a v . . 


¥ = hn'vnj- V v 2.i> + B i. ix ~ri i+ hi 


2.1 


'^ x ax 8 


flu 

'i3y a 


( 8 . 2 ) 


and the corresponding linearized system as 


5 v. . 2 

— -iiL - 2 b . , v, . + D . . 

at ic=1 ik kD idx 


a 8 y. 


y + Ih 


-2 
a v. 


ax 




ay 


¥ 


(8.3) 


where 


df. . 

b ik = V u s • 


The system (8,2) and (8,3) are associated with the 
following initial, boundary and matching conditions, 

i Initial conditions : 


v i;j (s>y,o) = P i3 (K) 


such that 


(8.4-a) 

(8,4-irb ) 
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(ii) Boundary conditions : 


v il 

(0,y 

II 

4^ 

(K 

0 

v iN 

(i >y ; 

,t) = 

0 

v. . 

(x,0, 

it) = 

0 

ID 


v. . 
ID 

(x,B, 

t) = 

0 


(iii) Matching conditions : 


V iD 



>t) 

(8.6 ) 

D n 

9v. - 

(-3^ ) = 

d. , . 

at x = L . - 
D- 1 

(8.?) 

for 

1=1,2 | 

j — 2 j « « • jN • 




In the following, we discuss the linear and nonlinear 
stability of equilibrium state in the cases of competition and 
prey predator models separately using the Liapunov second method. 

8.3 COMPETITION MODEL 


The interaction function for competing species including 
crowding effects can be written as follows: 


f l3 = U l3 (a l- b l U ir C l U 23 } 
f 2j = U 2^ a 2" b 2 U 23“ C 2 U lD^ 


( 8 . 8 ) 


which give the positive uniform equilibrium state for the 
systems (8.2), (8.3) as 


_ ^'l b 2"" a l a 2 

= b-, b 0 -c.c 0 

1 (>* X 


^2 ^ 1*"*^ 2^1 

U 2 ~ IdTTdTTcTcT 
12 12 
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provided 


zi > i > c i 
c 2 a 2 F“ 

(8.9) 

£± < fl < C 1 

C 2 a 2 ^ 

(8.10) 

°° rreS ^ nai “« linearised system is given by ( 8 .3) wlth 

b ll = “ b l a l < 0 » b 12 = “ C 1 S 1 < o 
b 2l = ~’ C 2 U 2 < 0 > h 22 = " b 2 S 2 < 0 

(8.11) 

8.3.1 DISCUSS IOIT Of LINEAR SIABILIIY 


lo discuss the linear stability, we choose 

the positive 

definite function e, as 

1 B 


E = 2 f f C v i + cv|3 dxdy 

0 o ^ 

(8.12) 

W hO 2? 0 V . rs tt ■pov? T - 

1 v i0 f0r L j-1 < X < L. t 

and 0 is a positive constant to be determined. 


(8.3) the time derivative of k is obtained as 

dE r B __ 

4t ' o jfl 1. _ x Lb U V lj + ° b S2* h +(b 12 +Cb Zl> 

V. .V 

10 22 

dv. . 2 0v o . 2 



(8.13) 


av.. 2 0v 2 

after using (8.11) and conditions (8.5)-(8.7). 
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It can be seen from (8.13) that in absence of dispersion, 
dE/dt is negative definite provided, 

b lV c l c 2 > 0 (8.14) 

obtained after using Sylvester conditions and choosing c=c ^/ c > 
The condition (8.14) assures the existence of Liapunov function 
in absence of dispersion and gives the necessary and sufficient 
conditions for stability in absence of dispersion. 

It can be noted from (8.13) that the equilibrium state 
which is stable otherwise remains so with dispersion in patchy 
habitat also. 

In case, when condition (8.14) i s not satisfied, let us 
explore the possibility for which E is still a Liapunov function 
with dispersion. 

Eor this, using the inequality (A10) obtained in the 
Appendix, the expression (8.13) for dE/dt simplifies to 

B L _ _2 2 _2 „ 

dt - £ £ B b n~ x l^ v l + °( b 22“ *2 ) v 2 + ^ b 12 +cb 2l^ V l V 2^ dxdy 

(8.15) 

which shows that dE/dt is negative definite provided 

_2 _2 _2 _2 

b ll b 22" b 12 b 2l" b ll X 2~ b 22 X 1 + X 1 X 2 > 0 (8.16) 

where 3/ ; i=l,2 are the roots of the equation obtained 
after simplifying (A6 ) and (A9 ) • It is noted that if the 

s 

2 

condition (8.16) is satisfied for minimum roots of • i=l p 

i 9 9 ~ 9 
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it would also be satisfied for all t no other roots, fhus for 

-2 

minimum values of X ^ i=l,2 (8.16) is the necessary and 

sufficient condition for linear stabilitv of otherwise unstable 
equilibrium state . 

8.3,2 DISOUSSIOM OF NONLINEAR STABILITY 


fo investigate the nonlinear stability of the equilibrium 
state (u^iig) of the system (8.2), the following positive 
definite function, in the positive quadrant of u^,Ug is 
considered , 


L v v n 

V = / / [v -u ln(l+ -±) + c { v -u ln(l+ ~)>)]dxdy 

oo ii u ^ £L 

(8.17) 

Using (8.11) we get the time derivative of Y as 



+ 


“iL-tx 

(S 1 + T 1;j ) 2 


av. . 2 u-jD-. dv.. 

( _H) + liL ( - M.) 
8x 


■h 


cu 2' D 3.ix 

( V v 2 j ) 


( !Li) 2 + °LLs_ 

ax 


3v 9 . 2 

(-af 2 ) H<^ay 


(8.18) 

after using boundary conditions (8.5) and matching conditions 
(8.6). It is seen from (8.18) that as in the case of linear 
stability, the condition for nonlinear stability of the 
equilibrium state is given by (8.14) in absence of dispersion. 
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Thus, the equilibrium state is nonlinearly stable if it is 
linearly stable without dispersion. Further it can also be 
noted from (8.18) that the otherwise stable equilibrium state 


remains stable nonline arly with dispersion also as 
are positive for i=l,2$ 3=1,2... I. 




Even if condition (8.14) is not satisfied, in the 
following, we find a subregion A of positive quadrant in 
which the equilibrium state (u^,ilg) becomes nonlinearly 
stable with dispersion. 


For this, let us assume 




3 = 1 , 2 ,...,! 


such that 


U l3 = V V 13 < U 03 

(8.19) 

U 2j = U 2 +V 2j < Y 03 


and U o . > u 1 , V 0 . > u 2 . 

Considering the region A = {(u 1 ,Ug) : u^ < U 0 , U g < > 

where U ,? are the maximum values of U . ,Y .j j=l,2,...! 
o o 03 03 

respectively . 


The expression (8.18)for dV/dt can now be simplified 
in the region A to the following form 
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Using the inequality (A10) we get from (8.20) 



+ (°l+ c c 2 ) y 2 ]dxdy 

which gives the following condition for the function 
he a Liapunov function 


b iV c i°2 + 


b l x 2 h 


b 2 X 1 1 
L 2 


+ 


^2 u ^ u g 

— 

Y o 


> 0 


( 8 . 21 ) 
Y to 


(8.22) 


where X^ and Xg are the roots of the equation obtained 

after simplifying (A6 ) by using matching conditions (A4) and 

boundary condition (A5). Thus, the sufficient condition for 

nonlinear stability of the equilibrium state in the region a 

2 

is given by (8.22). The equations for determining i=i,2 

are given by (All), (A12) in cases of habitats with two and 
three patches respectively. As discussed in Appendix, the 
minimum root of these equations increase *vith dispersion in 
patchy habitat, it can he concluded that the region of stability 
obtained from the condition (8.22) increases with increasing 
dispersion coefficients in patchy habitat. 

In a particular case when crowding terms in interaction 
functions for the two species are zero, i.-e. b^=bg=0, it can 
be noted from (8.11) and (8.14), that the equilibrium state 
which is unstable without dispersion may become linearly stable 
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provided 

_2 _ 2 , . 
X 1 Xg -o ± c 2 > 0 ( 8.23 ) 

[see equation (8.16)j. However, from (6.22), the region of 
nonlinear stability can now be obtained from the following 
inequality 

*1 *2 “ c l°2 > ° (8.24) 


Comparing (8.23) and (8.24) with (8.16), (8,22) it can be 
remarked that the effects of crowding with dispersion in a 
patchy habitat are to increase the degree of linear stability 
and region of nonlinear stability. 


Further it can be noted from (A9) that in the case of 
one dimensional habitat where D^y=0; i=l,2; 3 = 1 , 2 ,. ..H, the 
region of nonlinear stability is smaller than that in the two 
dimensional habitat. Ihus the equilibrium state which is 


nonlinearly stable in a two dimensional habitat may not be so 
in a linear one dimensional habitat . 


8.4 PBEX-PREDAIOR MODEL 

In case of prey predator model, the Vol terra type of 
interaction functions are given by 

f ij (u ir u 2D } = (a r b i u ir c i u 2D } Ul 3 

f 2o (u l d ,U 2D^ = (a 2" D 2 U 2r C 2 U l0 ) U 2D 

which gives the positive equilibrium state as 
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a 2 C l + 

a l b 2 

a l C 2~ 

a c o i 

£ x 

°iV 

b l b 2 ’ 

c l c 2+ 

b i b 2 


provided a^c^-a^b^ ^ -^ e linear lzeu system can 

be given 

by (8.3 ) with 



b ll = " b l U l ?b 12 = ‘ 

- c l^l 




(8.26) 

b 22 = ~ b 2 U 2 ’ b 2 1 = 

C 2 U 2 


To see the linear stability of 

the equilibrium 

state , 


the same positive definite function E as (8.12) is 
considered, from which, proceeding as before, it can be noted 
that E is a Liapunov function as the condition (8.14) in 
this case is automatically satisfied. Thus the equilibrium 
state which is linearly stable otherwise would remain so with 
dispersion in patchy habitat. She same result is true for 
nonlinear stability also and this can be verified by choosing 
same Liapunov function Y as (8.17). 

8.5 CONCLUSIONS 

In this chapter, the effects of patchines arising out of 
variable dispersion coefficients on linear and nonlinear 
stability have been studied by using Liapunov second method. 

In the case of competition model it has been shown that 
the otherwise stable equilibrium state remains linearly as 
well as nonlinearly stable with dispersion in patchy habitat. 
However, if the equilibrium state is unstable without dispersion 



177 


the necessary and sufficient conditions for the linear 
stability of the equilibrium state in the positive quadrant 
has been obtained. But in the case of nonlinear stability 
with dispersion it can be seen that the equilibrium state 
could only be stable in a subregion of positive quadrant 
and this region would increase with increasing dispersion. 

In case of prey-predator model, the otherwise stable equili- 
brium state remains stable linearly and nonlinearly with 
dispersion. 
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-APPENDIX 


Let us consider 


L 


L • 
3 


1 " A D i.ix { v li J dx > 0 


131 


id 


(A-l) 


where In.., i=l,2j j=l,2,...N are arbitrary differentiable 


functions in L < x < L . . 

D -1 “ D 


Ihe integral (Al) can be simplified to get 


N L j av. 4 2 


-L D ijx / 3x^ “ ^ h ij“ h i^ v i;j3 dx 


D-l 


D-l 


N 


+ ^J* l> lj (L 3 >f 3 (L 3 )-h i3 (L d . 1 )vf j (l 3 . 1 )J >C 


(A-2 ) 


Let h. . (x) to be the solution of the following differential 
*■ D 


equation 


W 3 (xJ - h8.(xj = ^ 


(A-3) 


where h! . is the derivative of h. . w.r.t. x and . being 
ID ID 3-D 


arbitrary real constant, subject to the conditions 


W« ( V = * 3=1.2, 

It = o 

x- o 

Lt T v iN^ x ^ b iN^ x ^ = 0 


. .H-l 

(A-4) 


(A-5) . 
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Then can be determined as follows: 

' h il = - x il cot x u 


^3 “ " X ij taI1 {X ij x + A ij } 


h iM 


- cot U, T<r (x-L) > 


iN 


U6) 


Now substituting (A3)-(A5) in (A2) we get the inequality. 
N L j av. .2 * r L 


N j 3v- .2 l\f o i o 

i=i D «x / <-#> ^ X h ( 

a- 1 J j-i 


dx 


2-1 

i=l,2 ; j = l,2, ...N 


(A-7) 


Thus using (A-7 ) , 
B N 


n im ^ 3 av. .2 BN L ;j a 2 v. . 

J 2 D. . J dxdy + f Z D. . J ( |^.)dxdy 


o ijj' ax 

B M 


o jir 1 ^ i . i ay 2 

J 

-r* 2 L . 

D, .„n _ 2 


> f Z [3 (D. . x.. + ) / vf. dx^jdy 

o i=l 13x 3-3 B* L. , 13 

0 i-l 


Assuming 


D. . X?. + hil = Xi l-WSi 3=1,2, ...H 

13X 13 b 2 


(A-8) 


(A-9) 


we get from (A8) 

B N L j 

/ S / C3>n 


3=1 L 


av, , 2 
(-TT 11 ) + 


i-1 


ijx v 3x 


av. . 2 

D 1 > 

iiT ay ; J - 


x i 


2 B L 

/ / 
o o 


v? dxdy 


(A-10) 


where can be determined from (A6) and (A9) using the 

conditions (A4). In the particular cases of habitats with 
two and three patches, the form of equations determining 
x. % i=l,2 can be written as 
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i - ll -j 

i 1 


cos C-— :IJ sinE 


VW 






n ( W-, 


Jsin 


x i JJ l- 


^x 


iOsC 


¥i -, 


JoosL 


!— 


:-TD lgx stoC-^jsmt;- 


^12x 




(JU12) 


respectively for equal dispersion coefficients in y-directions. 
'fhe above two equations have been solved numerically in 
Chapter IV and it has. been pointed out that minimum roots of 
these two equations increase with increasing dispersion 
coefficients in different patches. 
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